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The extremely high precision of current radio/optical interferometric observations and the unpar- 
alleled sensitivity of existing (LIGO) and future (LISA, ASTROD) gravitational-wave detectors 
demand a much better theoretical treatment of relativistic effects in the propagation of electromag- 
netic signals through variable gravitational fields. Especially important for future gravitational-wave 
observatories is the problem of propagation of light rays in the field of multipolar gravitational waves 
emitted by a localized source of gravitational radiation. A consistent approach giving a complete 
and exhaustive solution to this problem in the first post-Minkowskian approximation of General 
Relativity is presented in this paper. We derive a set of equations describing propagation of an 
electromagnetic wave in the retarded gravitational field of a time-dependent localized source emit- 
ting gravitational waves with arbitrary multipolarity and show for the first time that they can be 
integrated analytically in closed form. We also prove that the leading terms in observable relativis- 
tic effects depend exclusively on the values of the multipole moments of the isolated system and 
its time derivatives taken at the retarded instant of time on the null cone and do not depend on 
their integrated values. By making use of our integration technique we reproduce the known results 
of integration of equations of light rays both in a stationary field of a gravitational lens and in 
that of a plane gravitational wave, thereby establishing a relationship between our formalism and 
the approximations used by previous researches. The gauge freedom of our formalism is carefully 
studied and all gauge-dependent terms in the expressions for observable quantities are singled out 
and used for physically meaningful interpretation of observations. Two limiting cases of small and 
large values of the light-ray impact parameter, d, are elaborated in more detail. We explicitly show 
that in the case of small impact parameter the leading order terms for any effect of light propaga- 
tion in the field of an arbitrary multipole depend neither on its radiative nor on its intermediate 
zone contributions. The main effect rather comes from the near zone terms. This property makes 
much more difficult any direct detection of gravitational waves by astronomical techniques if general 
relativity is correct. We also present an analytical treatment of time delay and light-ray bending 
in large impact parameter case corresponding to the approximation of a plane gravitational wave 
of arbitrary multipolarity. Explicit expressions for time delay and deflection angle are obtained in 
terms of the transverse-traceless (TT) part of the space-space components of the metric tensor. 
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1. NOTATIONS AND CONVENTIONS 



Metric tensor on the space-time manifold is denoted by gap and its perturbation hap = Sap — "Hap- The determinant 
of the metric tensor is negative and is denoted as g — dct\\gap\\- The four-dimensional fully antisymmetric Lcvi-Civita 
symbol Sap-yS is defined in accordance with the convention £0123 = +1- 

In the present paper we use a geometrodynamic system of units [1-4] such that the fundamental speed, c, and 
the universal gravitational constant, G, are equal to unity, that is c = G = 1. Space-time is assumed to be globally 
asymptotically-flat and covered by a single coordinate system a;^, x^, a;'') = {t,x,y,z) , where t and {x,y,z) are 
time and space coordinates respectively. This coordinate system is reduced at infinity to the Minkowskian coordinates 
defined up to a global Lorentz-Poincare transformation [5]. We shall also use the spherical coordinates (r, 9, (f) related 
to (x, y, z) by the standard transformation 

X = r sin 6 cos , y = r sin 6 sin , z ~ r cos 9 . ( 1 ■ 1 ) 

Greek (spacetime) indices range from to 3, and Latin (space) indices run from 1 to 3. If not specifically stated the 
opposite, the indices are raised and lowered by means of the Minkowski metric rjap = diag(— 1, 1, 1, 1). Regarding this 
rule the following conventions for coordinates hold: = Xi and ~ —a'o- We also adopt notations, Sij = diag(l, 1, 1), 
for the Kroneker symbol (a unit matrix), and, e^jfe, for the fully antisymmetric 3-dimcnsional symbol of Levi-Civita 
with convention £123 = -1-1. 

Repeated indices are summed over in accordance with the Einstein's rule [1-3] . In the linearized approximation of 
general relativity used in the present paper there is no difference between spatial vectors and co-vectors as well as 
between upper and lower space indices. Therefore, for a dot product of two space vectors we have 

A'B, = A,B, = AiBi + A2B2 + . (1.2) 

In what follows, we shall commonly use spatial multi-index notations for Cartesian three-dimensional tensors [6], that 
is 



IAi — 1ai...ai- (1-3) 

Tensor product of / identical spatial vectors fc' will be denoted as a three-dimensional tensor having I indices 

kax...^ai — ^ai...ai ■ (-^■^) 

Full symmetrization with respect to a group of spatial indices of a Cartesian tensor will be distinguished with round 
brackets put around the indices 

Q{ai...ai) = J^'^Qa{l)...a(l) , (1-5) 

where a is the set of all permutations of (1, 2, /) which makes Qai...ai fully symmetrical in ai . . . a/. 

It is convenient to introduce a special notation for symmetric trace-free (STF) Cartesian tensors by making use of 
angular brackets around STF indices. The explicit expression of the STF part of a tensor Qai...ai is [6] 

[1/2] 

Q<ai...ai> — y ] 0'k^{aia2 ' ' ' ^a2k-ia2k^a2k+i---o.i)bibi...bkbk t (1-6) 
fe=0 

where [1/2] is the integer part of 1/2, 

Sai...a, = Q{ai...a,) (1-7) 

and numerical coefficients 

, _ (-1)*^ U {2l~2k-l)\\ 



^ (2fc)!! (2? - 1)!! {l-2k)\ 
We also assume that for any integer I > 



(1.8) 



n = l{l ~1)...2 -1 , 0!ee1, (1.9) 
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and 



One has, for example, 



= 2)(/- 4). ..(2 orl) , 0!! = 1 . (1.10) 



T<cabc> — Ti^abc) ^ '^^abT(cjj) ~ -^^bcT(ajj) ~ -^5acT{bjj) ■ (l-H) 

Cartesian tensors of the mass-type (mass) muhipolcs X<^,> and spin-type (spin) multipoles 5<^,> entirely describing 
gravitational field outside of an isolated astronomical system are always STF objects that can be checked by inspection 
of the definition following from the multipolar decomposition of the metric tensor perturbation hap [6] . For this reason, 
to avoid appearance of too complicated index notations we shall omit in the following text the angular brackets around 
indices of these (and only these) tensors, that is we adopt: Iai =^<Ai> and 5^1, = S^Ai>- 

We shall also use transverse (T) and transverse-traceless (TT) Cartesian tensors in our calculations [2, 6]. These 
objects are defined by making use of the operator of projection Pjk = 5jk — kjk onto the plane orthogonal to a unit 
vector kj. Thus, one has [6, 7] 

Qai...ai = PaibiPa2b2---PaibiQbi...bi 7 (1-12) 

[1/2] 

^l-^(aia2 ■ ■ ■ Pa.2k-ia:2k^a2k + i---ai)bibi...bkbk : (1-13) 



'ai...ai 

fc=0 



where again [1/2] is the integer part of 1/2, 

' {ax...ai) 



Wa,...a, ^ Ql^,,,,. (1.14) 



and numerical coefficients 



For instance, 



, (-l)M(/-fc-l)!! 

4^ k\il-2k)\ • ^^-^^^ 



Qll = \ {Pa^PbJQ^J + Pb^PaJQ^J) ~ ^Pab {PjkQjk) ■ (1.16) 

Polynomial coefficients will be used in some of our equations and they are defined by 

Ci{pi,...,p„) = —^ -, (1.17) 

PV ■ --Pn' 

where / and pi are positive integers such that J2"=iPi — ^- We introduce a Heaviside unit step function H{p~ q) such 
that on the set of whole numbers 

HiP-q)^\'; ''/^'^ (1.18) 

[1, ifp>q. 

For any differentiable function / = f{t,x) one uses notations: /q = df/dt and f i ~ df/dx' for its partial 
derivatives. In general, comma standing after the function denotes a partial derivative with respect to a corresponding 
coordinate: /,„ = d f (x) / dx°' . Overdot denotes a total derivative with respect to time / = df/dt = df/dt + x'^df/dx'^ 
usually taken in this paper along the light ray trajectory x{t). Sometimes the partial derivatives with respect to 
space coordinate x' will be also denoted as di = d/dx"^, and the partial time derivative will be denoted as dt = d/dt. 
Covariant derivative with respect to the coordinate will be denoted as Vq. 

We shall use special notations for integrals with respect to time and for those taken along the light-ray trajectory. 
Specifically, the time integrals from a function F{t,x), where x is not specified, are denoted as 

t t 
F^-^'>{t,x)= J F{T,x)dT , F^~^'>{t,x)= J F^-^'>{T,x)dT . (1.19) 



6 



Time integrals from the function F(t,x) taken along the light ray so that spatial coordinate x is a. function of time 
X = x{t), are denoted as 

t t 
F^-^\t,r) = J F{T,r{T))dT, F^-^\t,r)= J F^-^\T,r{T))dT . (1.20) 

— oo — oo 

Integrals in equations (1.19) represent functions of time and space coordinates. Integrals in equations (1-20) are 
defined on the light-ray trajectory and are functions of time only. 
Multiple time derivative from fimction F{t, x) is denoted by 

F"Ht.^)^'-^. (L21, 

so that its action on the time integrals eliminates integration in the sense that 

pi.) u^. d^^^Fi-^\t^x) dr'+^Fi^-Ht^x) 

and 

dP+^Fl-^\t,x) dP+^Fl-^\t,x) 

Spatial vectors will be denoted by bold italic letters, for instance A* = A, k'^ = k, etc. The Euclidean dot product 
between two spatial vectors a and b is denoted with dot between them: a^bi ^ a ■ b. The Euclidean wedge (cross) 
product between two vectors is denoted with symbol x, that is eijka^h^ = (a x b)'. Other particular notations will 
be introduced as they appear in the text. 



2. INTRODUCTION 



2.1. Statement of the Problem 



Direct experimental detection of gravitational waves is a fascinating, yet unsolved problem of modern fundamental 
physics [8-12]. Enormous efforts are applied to make a progress in its solution both by theorists and experimentalists. 
The main theoretical efforts are presently focused on calculation of templates of the gravitational waves emitted 
by coalescing binary systems comprised of neutron stars and/or black holes [13-18] as well as creation of improved 
filtering technique for gravitational wave detectors [19, 20] which will enable one to extract the gravitational wave 
signal from all kind of interferences present in the data collected by the gravitational wave observatories. Direct 
experimental efforts have led to the construction of several ground-based optical interferometers [21-25] with arms 
reaching few miles. Certain work is under way to build super-sensitive cryogenic-bar gravitational-wave detectors of 
Weber's type [26-29]. Laser interferometric space projects such as LISA [30] and ASTROD [31] aimed to significantly 
increase the sensitivity of the gravitational-wave detectors are currently under intensive discussion at NASA, ESA 
and China. 

No doubt, detection of gravitational waves by the specialized gravitational wave antennas would give the most direct 
evidence of the existence of these as yet elusive ripples of the space-time. At the same time it is conceivable that 
there exist a number of astronomical events which observation with a significantly improved astronomical technique 
might be used for indirect detection of gravitational waves as well as the extension of our knowledge of the conditions 
existing in a very early universe and of physics of the astrophysical systems emitting gravitational waves. It is worth 
emphasizing that the specialized ground-based and spaceborn gravitational wave detectors are sensitive to a limited 
band of gravitational wave frequencies ranging from 10^ Hz down to 1 Hz [31-33]. Hence, they can not tell us anything 
about the astrophysical effects produced by time-dependent gravitational fields having no radiative character but 
still playing an important role in the near and/or intermediate zones of isolated gravitating systems. Furthermore, 
ultra-long gravitational waves with frequencies much less than 1 Hz can be generated by localized gravitational 
sources in the early universe [34-38] and by ensembles of binary systems in our and other galaxies thus forming 
a statistical background with a reasonably large energy density [39-41]. Near-zone time-dependent perturbations 
of the gravitational field and the ultralong gravitational waves can be studied only via astronomical observations 
(for example, pulsar timing, radiometry of cosmic microwave background radiation, etc.) which provide us with an 
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additional window for exploring different properties of the gravitational field thereby extending the frequency range 
of the gravitational wave astronomy to all plausible frequencies [32, 33]. 

Any astronomical observation is done with electromagnetic waves (photons) of different frequencies with an overall 
spectrum spreading out from radio waves to gamma rays. Therefore, in order to detect the effects of the gravitational 
field by astronomical technique one has to elaborate on the problem of propagation of electromagnetic waves through 
this field generated by a localized gravitationally-bounded distribution of masses moving with velocities and acceler- 
ations determined from the equations of motion of a particular theory of gravity [42]. Equations of propagation of 
electromagnetic signals must be derived in the framework of the same theory of gravity in order to retain description of 
the gravitational and electromagnetic phenomena on a self-compatible ground [43] . In the present paper we rely upon 
the Einstein theory of general relativity (GR) and assume that light propagates in vacuum, that is the interstellar 
medium has no impact on the light propagation. 

General Relativity is a geometrized theory of gravity and it assumes that both gravity and electromagnetic field 
propagate locally in vacuum with the same speed which is equal to the fundamental speed of special theory of relativity 
[1-3]. Gravitational field is found as a solution of the Einstein field equations. The electromagnetic field is obtained 
from the Maxwell equations on a curved space-time derived in accordance with the principle of equivalence [1-3] which 
determine how the affine connection FJ^^ is related to the metric tensor g^^ . In the approximation of geometric optics 
[1-3] electromagnetic signals propagate along null geodesies of the metric tensor [1-3] . 

The problem of finding solutions of the null ray equations in curved space-time attracted many researchers since 
the time Einstein discovered the field equations of General Relativity. Exact solution of this problem was found in 
many, particularly simple cases of space-time possessing specific symmetries like the Schwarzschild or Kerr black holes, 
homogeneous and isotropic Friedman-Robertson- Walker cosmological model, plane gravitational wave, etc. [1-3, 44]. 
These exact description of null geodesies was (and in some cases still is) satisfactory for purposes of experimental 
gravitational physics for quite a long time. However, real physical space-time has no symmetries and solution of 
such current problems as direct detection of gravitational waves, interpretation of an anisotropy and polarization of 
cosmic microwave background radiation (CMBR), exploration of inflationary models of the early universe, finding 
new experimental evidences in support of general relativity and quantum gravity, development of higher-precision 
relativistic algorithms for space missions, and many others, can not fully rely upon the mathematical techniques 
developed exceptionally for dealing with symmetric space-times. Especially important is to find out a method of 
integration of equations of light geodesies in space-times having time-dependent gravitational perturbations. 

In the present paper we consider a general case of an asymptotically-flat space-time to which an isolated astronomical 
system is embedded. We assume that gravitational field of this system has a whole set of gravitational multipoles 
emitting gravitational waves [45] . Precise mathematical definition of the asymptotic flatness of space-time based on 
the concept of conformal infinity was worked out in [46-51] (see [3] for a comprehensive mathematical introduction 
to this subject). This definition implies conformal transformation technique and Hansen's multipole moments [52] 
which work perfectly in stationary space-times but are not applicable for analysis of real astronomical observations 
and for interpretation of relativistic effects in propagation of electromagnetic signals in time-dependent gravitational 
fields. We are going to study linear effects of gravitational multipoles on light propagation so that relativistic effects 
produced by non-linearity of the gravitational field will be ignored. Linearized approximation of general theory of 
relativity represents a straightforward approach for description of the multipolar structure of the gravitational field 
of a localized astronomical system. Multipolar formalism was developed by Thorne [6] and further extended and 
elucidated in [53-57]. 

Exact formulation of the problem under discussion in this paper is as follows (see Figure 1 for more detail). The 
space-time is covered by a single coordinate system = {ct,x), where t is time and x are space coordinates. An 
electromagnetic wave (photon) is emitted by a source of light at time to from a point Xq towards observer which 
receives the photon at time t at the point x. Photon propagates through the time-dependent gravitational field of 
an isolated astronomical system emitting gravitational waves. The structure of the gravitational field is described by 
a set of multipole moments of the system which are functions of the retarded time t ~ r/c where r is the distance 
from the isolated system to the field point and c is the fundamental speed. The retardation is due to the finite speed 
of gravity which coincides in general relativity with the speed of light in vacuum [58] . Our problem is to find out 
relationships connecting various physical parameters (direction of propagation, frequency, polarization, intensity, etc.) 
of the electromagnetic signal at the point of emission with those measured by the observer. Gravitational field affects 
propagation of the electromagnetic signal and changes its parameters along the light ray. Observing these changes 
allow us to study properties of the gravitational field of the astronomical system and to detect gravitational waves 
with astronomical technique. The rest of the paper is devoted to the mathematical solution of this problem. 
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2.2. Historical Background 

Propagation of electromagnetic signals through static and/or stationary gravitational field is rather well-known 
subject. Presently almost any standard textbook on relativity describes solution of the null geodesic equations in 
the field of a spherically-symmetric and rotating body. This solution is practical since it is used, for example, for 
interpretation of gravitational lensing events taking place in our galaxy and/or in cosmology [59-61]. However, 
continually growing accuracy of astronomical observations demands much better treatment of secondary effects in the 
propagation of light produced by the perturbations of the gravitational field associated with higher-order gravitational 
multipoles. Time-dependent multipoles emit gravitational waves which perturb propagation of light. Calculation of 
these perturbations is important for developing correct strategy for searching gravitational waves by making use of 
astronomical technique. 

Among the most interesting sources of periodic gravitational waves with a well-predicted behavior are binary 
systems comprised of two stars orbiting each other around common barycentcr. Indirect evidence in support of 
existence of gravitational waves emitted by binary pulsars was given by Taylor with collaborators [62. 63]. However, 
direct observation of gravitational waves remains an unsolved problem of experimental gravitational physics as yet. 
The expected spectrum of gravitational waves extends from ~ 10^ Hz to IQ-^^ Hz [32, 33]. Within that range, the 
spectrum of periodic waves from known binary systems extends from about 10"'^ Hz - the frequency of gravitational 
radiation from a contact white-dwarf binary [39], through the 10"'' to 10~^ Hz - the range of radiation from the 
main-sequence binaries [40], to the 10"^ to 10~^ Hz - the frequencies emitted by binary supermassivc black holes 
presumably lurking in active galactic nuclei [41]. The dimensionless strain of these waves at the Earth, h, may be 
as great as 10~^^ at the highest frequencies, and as great as 3 x 10~^'^ at the lowest frequencies in the spectrum of 
gravitational waves [32, 33]. 

Sazhin [64, 65] first suggested detection of gravitational waves from a binary system using timing observations of a 
pulsar, the line of sight to which passes near the binary. He had shown that the integrated time delay for propagation of 
the electromagnetic pulse near the binary is proportional to 1 /(P where d is the impact parameter of the unperturbed 
trajectory of the signal. More recently, Sazhin & Saphonova [66] have made estimates of the probability of observation 
of this effect for pulsars in globular clusters and showed that the probability can be high, reaching 0.97. Wahlquist [67] 
proposed another approach to the detection of periodic gravitational waves based on Doppler tracking of spacecraft 
travelling in space. His approach is restricted by the plane gravitational wave approximation developed earlier by 
Estabrook & Wahlquist [68]. Tinto ([69], and references therein) made the most recent theoretical contribution in 
this area. The Doppler tracking technique is used in space missions, by seeking the characteristic triple signature, 
the presence of which would reveal the infiuence of a gravitational wave crossing the line of sight from spacecraft to 
observer [70]. 

Quite recently, Braginsky et al. [71, 72] have raised a question about using very-long baseline radio interferometer 
as a detector of stochastic gravitational waves. This idea has also been investigated by Kaiser & Jaffe [73] and, 
in particular, by Pyne et al. [74] and Gwinn et al. [75] who showed that the overall effect is proportional to the 
strain of the metric perturbation caused by the plane gravitational wave. They also set an observational limit on the 
energy density of ultra-long gravitational waves present in early universe. Montanari [76] studied perturbations of 
polarization of electromagnetic radiation propagating in the field of a plane gravitational wave and found that the 
effects are exceedingly small. 

Fakir ([77, 78], and references therein) has suggested to use astrometry to detect periodic variations in appar- 
ent angular separations of appropriate light sources, caused by gravitational waves emitted by isolated sources of 
gravitational radiation. He was not able to develop a self-consistent approach to tackle the problem with necessary 
completeness and rigor. For this reason his estimate of the effect is too optimistic. Another attempt to work out 
a more consistent approach to the calculation of the light deflection angle by the field of arbitrary source of grav- 
itational waves has been undertaken by Durrer [79]. However, the calculations have been done only for the plane 
wave approximation and the result obtained was extrapolated for the case of the localized source of gravitational 
waves without justification. For this reason the magnitude of the light defiection angle was overestimated. The same 
misinterpretation of the effect can be found in the paper by Labeyrie [80] who studied a photometric modulation of 
background sources by gravitational waves emitted by fast binary stars. Because of this, the expected detection of 
the gravitational waves from the observations of the radio source GPS QSO 2022+171 undertaken by Pogrcbcnko et 
al. [81] was not based on firm theoretical ground and did not lead to success. 

Damour & Esposito-Farese [82] have studied the defiection of light and integrated time delay caused by the time- 
dependent gravitational field generated by a localized material source lying on the null direction being close to the 
line of sight. They worked in a quadrupole approximation and explicitly calculated effects of the full, retarded 
gravitational field in the near, intermediate, and wave zones by making use of the Fourier-decomposition technique. 
Contrary to the claims of Fakir [77, 78] and Durrer [79] and in agreement with Sazhin's [64] calculations, they found 
that the deflections due to both the wave-zone gravitational wave and the intermediate-zone retarded fields vanish 
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exactly. The leading total time-dependent deflection is given only by the quasi-static, near-zone quadrupolar piece of 
the gravitational field. 

Damour and Esposito-Farese [82] analyzed propagation of light in the quadrupolar field of the isolated system 
emitting gravitational waves under simplifying condition that the impact parameter of the light ray is small with 
respect to the distances from both observer and the source of light to the isolated system. We have found [7, 83] 
another way to solve the problem of propagation of electromagnetic waves in the quadrupolar field of the gravitational 
waves emitted by the system without making any assumptions on mutual disposition of observer, source of light, and 
the system, thus, drastically overpowering the result of Damour's paper [84]. At the same time the paper [83] did 
not answer the question about the impact of other higher-order gravitational multipoles of the isolated system on the 
process of propagation of electromagnetic signals. This might be important if effective emission of an octupole and/or 
higher multipolar gravitational radiation is equal or even exceeds that of the quadrupole as it may be in case of, for 
example, highly asymmetric star collapse [85], nearly head-on collision of two stars, or break-up of a binary system 
caused by recoil of two black holes [86] . 

In the present paper we work out a systematic approach to the problem of propagation of light rays in the field of 
arbitrary gravitational multipole. While the most papers on light propagation consider both a light source and an 
observer as being located at infinity we do not need these assumptions. For this reason our approach is more general 
and applicable for any disposition of the source of light and observer with respect to the source of gravitational 
radiation. The integration technique which we use for finding solution of the equations of propagation of light rays 
was invented in [7, 83]. 

The metric tensor and coordinate systems involved in our calculations are described in section 3 along with gauge 
conditions imposed on the metric tensor. The equations of propagation of electromagnetic waves in the geometric 
optics approximation are discussed in section 4 and the method of their integration is given in section 5. 

Exact solution of the equations of light propagation and the form of relativistic perturbations of the light trajectory 
are obtained in section 6. Section 7 is devoted to derivation of basic observable relativistic effects - the integrated time 
delay, the deflection angle, the shift of frequency, and rotation of the plain of polarization. We discuss in sections 8 and 
9 two limiting cases of the relative configurations of the source of light, the observer, and the source of gravitational 
waves - the gravitational-lens configuration and the case of a plane gravitational wave. Appendix gives approximate 
expressions for the Christoffel symbols, for the metric tensor taken on the light-ray trajectory, and for the gauge 
functions. 



3. THE METRIC TENSOR, GAUGES, AND COORDINATES 

3.1. The Canonical Form of the Metric Tensor Perturbation 

We consider an isolated astronomical system emitting gravitational waves and assume that gravitational field is 
weak everywhere so that the metric tensor is expanded in a Taylor series with respect to the powers of gravitational 
constant G. The metric tensor is a linear combination of the Minkowski metric, 77a/3j and a small perturbation hap 

gap = Vap + Gh^p + 0{G^) , (3.1) 
g"^ = Vcp - Gh^^ + 0[G^) , (3.2) 

where h^p 1 and we use T]ap to rise and lower indices so that, for example, h"^ = rii°'^rj^"h^^. For computational 
convenience and simplification of physical interpretation of our results, we put the origin of the coordinate system to 
the center of mass of the isolated astronomical system under consideration [87]. 

The most general expression for the linearized perturbation of the metric tensor outside of the astronomical system 
under consideration was given by Thorne [6] (see also [53]) in terms of its symmetric and trace- free (STF) mass and 
spin multipole moments. This perturbation is described by the following expression 

haP = h^^p - + Wa,p + Wp^a , (3.3) 

where Wa are, the so-called, gauge functions. The canonical perturbation h'^^' obeys to vacuum wave equation 



(3.4) 
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which solution is chosen as 



i^can. 



i^can. 



2M ^{-ly \lA,{t-r) 
' r +^ II 

1=2 



2eipgSp{t-r)Ng _ {-lYl \ £ipqSpA,_, {t - r) 

{i + iy. 



4E 



(-1) 

n 



1=2 



(3.5) 
(3.6) 



,qAi 



7 can. 



1=2 

-X. ./,can. , can. 



-I ,Ai 



(-1)' 



1=2 



II 



1ijAi^2 ^ ^) 



(-1)'/ 



^pq{iSj)pAi_2i^ ^) 



(3.7) 
(3.8) 



,qAi. 



Here Ai and Si are total mass and spin of the system, and Xai and Sai are two independent sets of mass-type and 
spin-type multipole moments, N'^ = jr is a unit vector directed from the origin of the coordinate system to the field 
point. Because the origin of the coordinate system has been chosen at the center of mass, the expansions (3.5) - (3.8) 
do not depend on the mass-type linear multipole moment 1; which is equal to zero by definition. We emphasize that 
in the linearized approximation the total mass M. and spin Si of the astronomical system are considered as constant 
while all other multipoles are arbitrary functions of time. It is well established however [1-3] that gravitational waves 
emitted away by the system reduce its energy as well as its linear and angular momenta. Hence, in higher order 
approximations Xi ^ 0, and the mass, spin, and linear momentum of the system must be considered as functions of 
time like any other multipole. 

The canonical metric tensor (3.5)-(3.8) depends on the multipole moments ^Aiit — r) and SAiit — r) taken at the 
retarded instant of time. The retardation is explained by the finite speed of propagation of gravity. In the near zone 
of the isolated system the retardation due to the propagation of gravity is small and all functions of time in the metric 
tensor can be expanded in Taylor series around time t [6, 53]. This near-zone expansion of the metric tensor can 
be smoothly matched to the solution of the linearized Einstein equations in the domain of space being occupied by 
matter. The matching helps to express the multipole moments in terms of matter variables [54] 
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In the first post-Newtonian approximation the multipole moments have a matter-compact support [54] 
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In the higher post-Newtonian approximations the multipole moments have no matter-compact support and are ex- 
pressed by means of more complicated functionals [88] . Radiative approximation of the canonical metric tensor reveals 
that a contribution of tails of gravitational waves must be added to the definitions of the multipole moments (3.10), 
(3.11). In other words, the multipole moments in the radiative zone of the isolated system read [55, 56, 89] 
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where 6 is a normalization constant which value is supposed to be absorbed to the definition of origin of time scale in 
the radiative zone but this statement has not been checked so far. 
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3.2. The Harmonic Coordinates 



Equation (3.3) holds in an arbitrary gauge imposed on the metric tensor. The harmonic gauge is defined by the 
condition 

2/i"^^ -h'" = 0, (3.14) 

where h = h"^. The gauge condition (3.14) reduces the Einstein vacuum field equations to the wave equations on 
gravitational potentials h"^ (see equation (3.4)). Harmonic coordinates x" are defined as solutions of the homogeneous 
wave equation □ = up to the gauge functions w°' . In particular, the harmonic canonical coordinates are defined 
by the condition that all gauge functions Wa = 0. The canonical metric tensor (3.5) - (3.8) depends on only two sets 
of multipole moments [6, 53] which reflects existence of only two degrees of freedom of free (detached from matter) 
gravitational field in general relativity [1-3]. At the same time one can obtain a generic expression for the harmonic 
metric tensor by making use of infinitesimal coordinate transformation 

x'" = - (3.15) 

from the canonical harmonic coordinates to arbitrary harmonic coordinates x'" with the harmonic gauge functions 
which satisfy to a homogeneous wave equation 

□ w"=0. (3.16) 
The most general solution of this equation contains four sets of STF multipoles [6, 53] 
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where Wa, , 3^iAi_i, and ZqAi_i are Cartesian tensors depending on the retarded time. Their specific form is 
a matter of computational convenience for derivation and interpretation of observable effects but it does not affect 
invariant quantities like the phase of electromagnetic wave propagating in the field of the multipoles. 

The most convenient choice simplifying the form of the metric tensor perturbations, is given by the following gauge 
functions 
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These functions substituted to equation (3.3) transform the canonical metric tensor perturbation to another, re- 
markably simple form 

2M 



hot 



p can. 



(3.21) 

(3.22) 
(3.23) 
(3.24) 
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where the TT-projection differential operator Pijki, apphed to symmetric tensors depending on both time and spatial 
coordinates, is given by 



(3.25) 



and A and A ^ denote the Euclidean Laplacian and the inverse Laplacian respectively. 

When comparing the canonical metric tensor with that given by equations (3.21)-(3.24) it is instructive to keep in 
mind that '^~'^^Xa, ~ = (t — r) and A(Xyi, (t — r) /r) = Xa, (t — r)/r for r 7^ 0. This is a consequence of the fact 
that function ^^^^Xa, {t — r) is a solution of the homogeneous d'Lambert's equation, that is, □ \^~'^^Iai {t — r)/r\ = 
for r ^ 0. We also notice that Xa,(< - r)/r = A"! X^, (t - r)/r and ("2)2:^^ -r)/r = A~i [lA,it - r)/r]. The 
metric tensor perturbation (3.21)-(3.24) is similar to the Coulomb gauge in electrodynamics [90]. 



3.3. The Arnowitt-Deser-Misner Coordinates 



The Arnowitt-Deser-Misner (ADM) gauge conditions in the linear approximation read [50, 91] 

2/ioi,i - = , 3/iijj - hjj^i = , (3.26) 

where for any function / = f{t,x) one uses notations: /,o = df/dt and f i = df /dx^. For comparison, the harmonic 
gauge conditions (3.14) in the linear approximation, read: 

2/ioi,i — hiifi = hoofi , 2/i,yj — hjj^i = — /loo.i ■ (3.27) 
The ADM gauge conditions (3.26) brings the space-space component of metric to the form 

grj = S,,il + ^hkk) + hf;^ , (3.28) 

where hfj^ denotes the transverse-traceless part of hij and hkk = 3/ioo. The ADM and harmonic gauge conditions 
are not compatible inside the regions occupied by matter. However, outside of matter they can co-exist. Indeed, it 
is straightforward to check out that the metric tensor (3.21)-(3.24) satisfies both the harmonic and the ADM gauge 
conditions in the linear approximation and assumption that Ai = [92]. We call the coordinates in which the metric 
tensor is given by equations (3.21)-(3.24) as the ADM-harmonic coordinates. 

The experimental problem of detection of gravitational waves is reduced to the observation of motion of test 
particles in the field of the incident gravitational wave. These test particles are photons of the electromagnetic signals 
and observers (mirrors in ground-based gravitational-wave detectors) which receive these photons. Gravitational 
wave interacts with both photons and observers and perturbs their motion. These perturbations must be explicitly 
calculated and clearly separated to avoid possible misinterpretations of observable effects of the gravitational wave. It 
turns out that the canonical form of the metric tensor (3.5) - (3.8) is extremely well-adapted for performing analytic 
integration of equations of light rays. At the same time freely- falling observers experience influence of gravitational 
waves emitted by the isolated astronomical system and move with respect to the canonical harmonic coordinate 
system in a complicated way. For this reason, observable effects imposed by the gravitational waves on the light 
propagation get mixed up with the motion of observer in these coordinates. 

Arnowitt, Deser and Misner [50] showed that there exist a canonical ADM coordinate system which has a special 
property such that freely-falling observers are not moving with respect to this coordinates despite that they are 
perturbed by the gravitational waves. This means that the ADM coordinates themselves are not inertial and, although 
have nice mathematical properties, should be used with care in the interpretation of physical experiments. Making use 
of the canonical ADM coordinates simplifies analysis of the gravitational wave effects observed at gravitational wave 
observatories (LIGO, LISA) or by astronomical technique because the motion of observer (proof mass) is excluded 
from the equations. However, structure of the metric tensor in such canonical ADM coordinate system does not allow 
us to integrate equations for light rays analytically because it contains terms that arc instantaneous functions of time 
[83]. 

The ADM-harmonic coordinates have advantages of both harmonic and ADM coordinates. Thus, the ADM- 
harmonic coordinates allow us to get a full analytic solution of the light-ray equations and to eliminate the effects 
produced by the motion of observers with respect to the coordinate grid caused by influence of gravitational waves. 
In other words, all physical effects of the gravitational waves are contained in the equations of light propagation. 
This conclusion is, of course, valid in the linear approximation of general relativity and is not extended to the second 



13 



approximation where gravitational-wave effects on light and motion of observers can not be disentangled and have to 
be analyzed together. 

Similar ideology based on introduction of the TT coordinate system is applied for analysis of the output signal of 
the gravitational-wave detectors with freely-suspended masses [1, 2, 26] placed to the field of a plane gravitational 
wave, that is at the distance far away from the localized system emitting gravitational waves where the curvature of 
the gravitational wave front is negligible. Our ADM-harmonic coordinates are essential generalization of the standard 
TT coordinates because it can be used at arbitrary distances from the isolated system emitting gravitational waves 
including its near and intermediate zones. 



4. EQUATIONS OF PROPAGATION OF ELECTROMAGNETIC SIGNALS 
4.1. Maxwell Equations in Curved Space-Time 

In this section all indices are raised and lowered by means of the metric tensor gap with g°'^ defined in accordance 
with standard rule g'^^gp-y = S"^. The general formalism describing the behavior of electromagnetic radiation in 
an arbitrary gravitational field is well known [93]. An electromagnetic field is defined in terms of the (complex) 
electromagnetic tensor Fafs as a solution of the Maxwell equations. In the high-frequency limit one can approximate 
the electromagnetic tensor Fap as [1, 2] 

Fap = Re{Aai3 cxp{iip)} , (4.1) 

where A^p is a slowly varying (complex) amplitude and (ys is a rapidly varying phase of the electromagnetic wave 
which is considered as a real function. Both amplitude and phase are functions of time and space coordinates. 
The source-free Maxwell equations are given by [1-3] 

V„i^^^ + VisF^a + V^F„/3 = , (4.2) 

VpF^"^ = , (4.3) 

where Vq denotes covariant differentiation [94]. Taking a covariant divergence from equation (4.2), using equation 
(4.3) and applying the rule of commutation of two covariant derivatives from the tensor field of a second rank [1-3], 
we obtain the covariant wave equation for the electromagnetic field tensor 

ngFai3 + Rcpi^.F^'" + R,.aF/ - R,.pFa^ = , (4.4) 

where = V"Vq, RaiSfiv is the Riemann curvature tensor, and Rap = R'^^-yp ^^^^ Ricci tensor (definitions of the 
Riemann and Ricci tensors in this paper are the same as in the textbook [3]). We consider the case of the gravitational 
field in vacuum where Rap = 0. Hence, in our case the equation (4.4) is reduced to more simple form 

dgFaP + i?a/3M.F^" = . (4.5) 

Differential operator in equation (4.4) taken along with the Riemann and Ricci tensors is de Rahm's generalized 
topological d'Alembertian [2, 95] for the electromagnetic field. 



4.2. Maxwell Equations in the Geometric Optics Approximation 

Let us now assume that the electromagnetic tensor Fap shown in equation (4.1) can be expanded with respect to 
a small dimensionless perturbation parameter e = l/L where I is a characteristic wavelength of the electromagnetic 
wave and L is a characteristic radius of space-time curvature. More specifically, we assume that the expansion of the 
electromagnetic field given by equation (4.1) has the form [2] 

-Fq/3 = (aQ/3 +e&Q/3 +e^Ca/3 + •■•) exp • (4-6) 

Substituting the expansion (4.6) into equation (4.2), taking into account the definition of the electromagnetic wave 
vector la = dip/dx°', and arranging terms with similar powers of e lead to the chain of equations 



la a/37 + + ^ Oq/S — , 

Vq, ap^ + \7p a^a + V7 Gap — {la bpj + Ip bja + bap) , 



(4.7) 
(4.8) 
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where we have neglected the effects of curvature which are of order and thus can be considered as neghgibly small. 
Similarly, equation (4.3) gives a chain of equations 

Ipa"^'^ = 0, (4.9) 

Vpa''^ + ilpb''^ = 0, (4.10) 

where we again neglected the effects of curvature. 

Equation (4.9) implies that the amplitude of the electromagnetic field tensor is orthogonal in the four-dimensional 
sense to vector la at least in the first approximation. Contracting equation (4.7) with la and accounting for (4.9), we 
find that the wave vector la is null, that is 

Za/" = 0. (4.11) 

Taking the covariant derivative of this equality and using the fact that la] = since la = Vq(/3, one can show that 
the vector la obeys the null geodesic equation 

fVpr = , (4.12) 
which means that the null vector is parallel transported along itself. Equation (4.12) can be expressed as 

_ + r^^//3r = o, (4.13) 

where A is an afhnc parameter along the light-ray trajectory, and 

= ^5"^ {d^g^p + dpg^^ - d,gp^) , (4.14) 

are the Christoffel symbols. 

Finally, equation (4.4) (or (4.5)) can be used to show that in the first approximation 

Dxaap + ^aap^O , (4.15) 

where D\ = Z^V^, and 

79=(l/2)V«r, (4.16) 
is the expansion of the light-ray congruence defined at each point of space-time by the derivative of the wave vector 

Equation (4.15) represents the law of propagation for the electromagnetic field tensor along the light ray. In general 
case of 1? ^ the electromagnetic field tensor is not parallel-transported along the light ray. It can be shown that 
the expansion ■)? of the light-ray congruence is defined by the stationary components of the gravitational field of the 
isolated astronomical system determined by its mass M , and spin 5' , and does not depend on higher-order multipole 
moments. It means that gravitational waves do not contribute to the expansion of the light-ray congruence in the 
linearized approximation and their impact on i!} is postponed to the terms of the second order of magnitude with 
respect to the universal gravitational constant G. 

4.3. Electromagnetic Eikonal and Light-Ray Geodesies 

1. The Congruence of Light Rays 

We shall assume that geometric optics approximation is a valid approximation and electromagnetic signals propagate 
in vacuum. This means that we consider electromagnetic signals with wavelengths much smaller than characteristic 
wavelength of gravitational waves emitted by the isolated astronomical system and that the speed of light is equal 
to the speed of propagation of gravitational waves. We also neglect all relativistic effects associated with explicit 
appearance of the curvature tensor of space-time in equations of light propagation. In accordance with consideration 
given in the previous section 44.2, we are allowed to describe electromagnetic signals as massless particles (photons) 
moving along light-ray geodesies with the initial-boundary condition [96] 

x{to) = xo, ' = k . (4.17) 
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These conditions determine spatial position Xq of the electromagnetic signal at the time of its emission, to, and 
direction of its propagation given by a unit vector, k, at the past null infinity, that is at the infinite distance and at 
the infinite past [2, 3] where the space-time is assumed to be flat (see Figure 1). We imply that vector k is directed 
towards observer. 

Several functions characterize tensor of the electromagnetic field Fap. We start in the next section from discussing a 
differential equation for the phase if of electromagnetic wave, which is also known as eikonal [1, 97]. In the next sections 
we will derive equations for the null geodesies and the Stokes parameters characterizing polarization properties of the 
electromagnetic wave. Exact differential equations for these quantities are too complicated and can not be solved 
exactly. For this reason we shall appeal to the approximations retaining in those equations only those terms which 
are linear with respect to the metric tensor perturbation hap- 

At the first iteration one can neglect relativistic perturbation in photon's motion and approximates it by a straight 
line 

=xiv(i) =4 + fc'(i-io) , (4.18) 

where Iq is the time of emission of the photon, Xq are space coordinates of the source of the electromagnetic signal 
taken at the time to, and fc' = fc is the unit vector along the trajectory of photon's motion defined by equation (4.17). 

The bundle of light rays makes 2+1 split of space by projecting any point in space onto the plane being orthogonal 
to the bundle (see Figure 2). This allows to make a transformation to new independent variables t and ^* defined as 
follows 

T^hx'j, , C^P)x'r,, (4.19) 

where — 5j — k'-kj is the operator of projection on the plane orthogonal to fc. It is easy to see that the parameter 
T is equivalent to time 

T = k-x = t-t* , (4.20) 

where t* = k ■ Xq — to is the time of the closest approach of the electromagnetic signal to the coordinate origin which 
is taken in our case coinciding with the center of mass of the isolated astronomical system. Because for each light 
ray the time t* is fixed one concludes that the time differential dr = dt. The reader may think that results of our 
calculations will depend on parameters and t* . This is, however, not true since and t* depend on the choice of 
the origin of the coordinate system and, thus, are unphysical. Inspection of resulting equations shows that parameters 

and t* vanish from the observed quantities as it must be in agreement with the physical point of view. 

The unperturbed light-ray trajectory (4.18) written in terms of the new variables (4.19) reads 

xiv(r) = fcV + e, (4.21) 

so that the new variable = ^ = fc x (a; x fc) should be understood as a vector drawn from the origin of the coordinate 
system towards the point of the closest approach of the ray to the origin. As vectors fc' and ^' are orthogonal, the 
unperturbed distance — ^/xixf between the photon and the origin of the coordinate system 



rN 



= Vt^ + (P , (4.22) 

where d = \^\ is the impact parameter of the unperturbed light-ray trajectory with respect to the coordinate origin. 

We introduce two special operators of partial derivatives with respect to r and determined for any smooth 
function taken on the light ray. These operators are defined as 

SO that, for example, 

k'di = . (4.24) 

An important consequence of the projective structure of the bundle of light rays is that for any smooth function 
F{t,x) defined on the light-ray trajectory one has 

F(t, x) 



\(- 


+ '4) 


\dx^ 


-fd 


ox^ J 





Pit,x) 



x—XQ-\-k{t — tQ 



x—xo-\'k{t — tQ) 



x—XQ-\--k{t — tQ) 



^,+h^)F{t*+r,^ + kr), (4.25) 
:|:F(t*+r,^ + fcr), (4.26) 
=^P(r +T,| + fcr). (4.27) 
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Here, in the left sides of Eqs. (4.25)-(4.27) one must first calculate partial derivatives and only after that substitute 
the unperturbed trajectory x = Xq + k{t — to), while in the right side one first substitute the unperturbed trajectory 
parameterized by variables r and ^* and then differentiate, ft is worth noticing that making use of equation (4.27) 
allows us to re- write equation (4.25) as follows 



dF{t,x) 



x—xo-\-k{t — to) 



W^^'^r'^^'W*^ + r , 4 + fcr) . (4.28) 



This equation is used later for decomposition of STF spatial derivatives from the retarded Lienard-Wiechert potentials 
of the gravitational field. 

2. The Eikonal Equation 

The eikonal is related to the wave vector of the electromagnetic wave as la = daf- This definition along with 
equation (4.11) immediately gives us a differential equation for the eikonal 

We are interested in an unperturbed solution of this equation which is a plane electromagnetic wave 

•^N ~ Vo+ UjkaX" , (4.30) 

where ipQ is constant, tu = v^o is (constant) frequency of the electromagnetic wave at infinity and ka is an 

unperturbed direction of the co-vector la- Equation (4.29) assumes that ka is a null co-vector with respect to the 
Minkowski metric in the sense that 

ll'^^kakp = . (4.31) 

We postulate that the co- vector ka = (— l,fc) where the unit Euclidean vector fc is defined at past null infinity by 
equation (4.17). 

In the linearized approximation of general relativity the eikonal can be decomposed in a linear combination of 
unperturbed, (^at, and perturbed, -01 parts 

(p ~ + LUip , (4.32) 

so that the wave co- vector 

+ ' ^^-^^^ 

Making use of equations (3.2) and (4.29)-(4.33) yields differential equation for the perturbed part of the eikonal 

This equation can be solved in all space by the method of characteristics [98] which, in the case under consideration, 
are the unperturbed light-ray geodesies given by Eq. (4.21). Hence, after making use of relationship (4.26) one gets 
an ordinary differential equation for finding the eikonal perturbation 

^ = \hT^-{r,m''kP + dr {k^w^ - w') , (4.35) 

where the canonical metric tensor perturbation h'^ap '^^ taken on the light-ray trajectory as shown in equations 
(B.1)-(B.9) and the gauge functions can be taken arbitrary. Equation (4.35) can be solved analytically with the 
mathematical technique shown in section 5. 
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3. Light-Ray Geodesic Equations 

General form of the geodesic equation for light rays is obtained from equation (4.13) after substitution of the 
approximate expressions (A.1)~(A.6) for the Christoffel symbols. It gives 

^'(0 =^'^00.1 - hoifi - i/ioo,oi* - hikfix'' - {hoi^k - ^ofc,i)i''- (4.36) 



hoo.kx''x^ - [htk,j - i^hkj.i^ x^x^ + (^hkjfl - hak,j^ x''x^x\ 



where overdot denotes an ordinary time derivative and /iqo, ft-oii hij are components of the metric tensor taken on the 
unperturbed light-ray trajectory as shown in equations (B.1)-(B.9), that is haf3 = haf3{t,XN{t)). 

Equation (4.36) can be further simplified after substituting the unperturbed light-ray trajectory (4.21) to the right 
side of Eq. (4.36) and making use of equation (4.25). Working in arbitrary coordinates one obtains 

= ^k'^k" d.h'i^^- - dr (^k^KT- - \ k'k^kPq^;''-^ - drr{w' - Vw°). (4.37) 

where all functions in equation (4.37) are taken (before any differentiation) on the unperturbed light-ray trajectory 
given by equation (4.21) and the gauge functions it;" have not yet been specified [99] which means that equations 
(4.37) are gauge-invariant. 

The main advantage of Eq. (4.37) is that when one integrates along the light-ray path the following rules, applied 
to any smooth fimction F{t,^), can be used 

J ^F{t, I) dr =F{t, + C{0 , (4.38) 



/ ^^(^' ^) dr =^ I F{t, dr . (4.39) 

This means that terms which are represented as partial derivatives with respect to the parameter r can be immediately 
integrated by making use of (4.38). At the same time Eq. (4.39) shows that one can change the order of operations 
of integration and differentiation with respect to the parameter which definitely simplifies the problem of finding 
solution of Eq. (4.37) as it will be evident in next sections. 

The gauge functions w" are not the only gauge functions which appear in the equations of motion (4.37) of the light 
ray. In order to see it more evidently, it is worth noticing that Eq. (4.37) is linear with respect to the perturbations 
of the metric tensor hafj. Hence, it can be linearly decomposed in two equations for perturbations of the light-ray 
trajectory caused separately by mass and spin multipole moments. Substitution of the metric tensor (3.5) - (3.8) to 
Eq. (4.37) and replacement of spatial derivatives with respect to with those with respect to parameters and t 
by making use of (4.27) yield 

^ ~ ^(G) ^" ^{M) ~^ ^{S) I (4.40) 

where x'j^j and i'^j are components of the photon's coordinate acceleration caused by mass and spin multipoles of 
the metric tensor respectively, and ij^) is the gauge-dependent acceleration. These components read 

= drr [{k'^U - ^U) + {k^w° - W^)] , (4.41) 



= 2(9, - hdr)— + ^d,Y,Y.Y. ^'(' -P.P-I, Q)H{2 -q)x 

1=2 p=0 (;=0 

- ^—J-^) ( 1 - J—jjk<ai...aJap^^...ai>df,'''d^ 



(4.42) 
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and 



CO i— 1 p 



1=2 p=0 q=0 



(^ + 1)! 



Q_iG-p-l,p-g,q)if(2-g)x 



(4.43) 
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where overdot denotes a time derivative which is equivalent to a partial derivative with respect to the parameter t* 
[100], are gauge functions generated by arbitrary post-Minkowskian coordinate transformations, and ip" are the 
gauge functions generated by the multipolar structure of the metric tensor in the course of reshuffling terms in the 
right side of equations (4.37) for a light-ray geodesic. 

More specifically, functions tp" have been obtained after we singled out in the equations of light propagation all 
terms containing second and higher-order time derivatives with respect to the parameter r. These terms enter the 
equations of motion in the form of combination fe'iys" — (p^ given in Appendix C. It is worthwhile to emphasize that 
we do not intend to separate /c^t/j" — (/j* in two functions cp'' and (/j" because such a separation is not unique while the 
linear combination fc'iyj" — (p' does. We have not combined functions (p" with the gauge functions for two reasons: 

1. to indicate that solution of equations of light-ray geodesic, performed in one specific coordinate system, leads 
to generation of terms which can be eliminated by gauge transformation, 

2. to simplify the final form of the result of the integration as all terms with second and higher order time derivatives 
are immediately integrated in accordance with Eq. (4.38). 

Gauge functions w" are arbitrary making our equations gauge-invariant. However, for the sake of physical interpre- 
tation of the result of integration of equations of light-ray geodesic, we shall chose function w"' to make our coordinate 
system both harmonic and ADM. Coordinate description of motion of free-falling particles in these coordinates post- 
pones effects of perturbations by gravitational waves emitted by multipoles to higher orders of approximation. Specific 
form of the gauge functions at arbitrary field point is shown in equations (3.17), (3.18) and their form at any 
point on the light-ray trajectory is given in Appendix C. 

It is important to realize that for each multipole of order I the right sides of Eqs. (4.42), (4.43) have such structure 
that all terms, which are proportional iol — th time derivative d\» from this multipole, mutually cancel out. Specifically 
for this reason all time integrals, which can not be performed explicitly, vanish from the relativistic perturbations 
of the light-ray trajectory irrespectively of the order of the multipole moment. This fact was established in [83] in 
a quadrupole approximation only. Present paper proves that this result is valid for arbitrary multipole. As one will 
see later this property of the null geodesic prevents to amplify effects of gravitational waves acting on a light ray 
propagating near a gravitational lens emitting multipolar gravitational radiation (binary star, etc.). 



4.4. Equations for Polarization of Light and the Stokes Parameters 



1. Reference Tetrad Field 



Algebraic classification of electromagnetic fields and studying their physical properties are conveniently performed 
with the help of the Newman-Penrose formalism [101, 102]. This formalism is based on introducing at each point of 
space-time a null tetrad of four vectors associated with the bundle of light rays defined by the electromagnetic wave 
vector P. The Newman-Penrose tetrad consists of two real and two complex null vectors (P, n", m", to"), where the 
bar indicates complex conjugation. The tetrad vectors are normalized in such a way that n^Z" = — 1 and mcjff' = +1 
are the only nonvanishing products among the four tetrad vectors. 
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The vectors of the null tetrad are not uniquely determined by specifying Indeed, for a fixed direction /" the 
normalization conditions for tetrad's vectors are preserved under the linear transformations (null rotation) [101, 102] 

l"^ = , (4.44) 

n'" = A-^ (n" + B?7i" + Bm" + BBl") , (4.45) 

m'" = e-^®(m" + Sr) , (4.46) 

m'" = e*® (m" + B^) , (4.47) 

where A, Q are real and i? is a complex scalar. These transformations form a four-parameter subgroup of the Lorentz 
group. 

For analysis of intensity and polarization of electromagnetic waves it is useful to introduce a local orthonormal 
reference frame of observer moving with four velocity and seeing the electromagnetic wave travelling in the positive 
direction of z axis of the reference frame. It means that at each point of space-time the observer uses a tetrad frame 

~ ( "^"o)' ^"i)' ^"2)' ^"(3) ) defined in such a way that 



e — u , e 



i-Lu'^)-' [r + (/^u'5)w"J , (4.48) 



and two other tetrad vectors, e"-^^ and 6^2)1 are unit space-like vectors being orthogonal to each other as well as to 
e"g^ and e"^-^ [103]. In other words, vectors of the tetrad are subject to the following normalization conditions 

9c.pel,)e%^ = V,. , v'-'eiy^,^ = 5"^ • (4.49) 

Let us define at each space-time point a coordinate basis 

^("o) = (1,0,0,0) , = (0,a\a2,a=') , d^^^ ^ {0,b\b'^,b^) , d^-^j ^ {0, k\ , k^) . (4.50) 

where the spatial vectors a — (a^, a^, a^), 6 = {b^,b^,b^), and k = {k^, k'^,k^) are orthonormal in the Euclidean sense, 
and the vector k defines the spatial direction of propagation of the light ray at infinity (see equation (4.17)). In a 
particular case of a linearized gravity field and static observer, when gap = ^]a0 + hap and u'^ = 1 — (l/2)/ioo, = 0, 
the locally-orthonormal tetrad E^^p-j is given by 

^^"o) = (1 + ^'^oo ,0,0,0^ , (4.51) 
a) = f . a' - Ih.jaA , (4.52) 



' (1) - yoj" . - 2' 
1 

— J 

2 



^"2) = hojb^ - " 7;huV ) , (4.53) 



= [h,,k= , F - \h,,k^^ . (4.54) 

which is a direct consequence of equations (4.48)-(4.50). The local tetrad e"^^ of the moving observer relates to the 
tetrad of the static observer by means of the Lorentz transformation 

e^/j) = >^''pEl^) , Ey^ = ~Xy"^^^ , (4.55) 

where the matrix of the Lorentz transformation is [2] 

A"o = u" = j (4.56) 

A"o = a'"' = , (4.57) 



u u 



X) = S'^ + —— , (4.58) 
1 -I- 7 

and the inverse matrix of the Lorentz transformation A"^ is obtained from A"^ by replacing m' — w*, as it follows 
from its definition. 
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The connection between the null tetrad and the frame e"^^ is given by equations 

r = _(/^u'')(e^o)+e^3)) , (4.59) 



-i(^«'')(e^o)-e'^3)) , (4.60) 



1 / 

1 / 



(e^i)+*e^2)j , (4.61) 



- *e^2) j • (4.62) 



If vectors e"^^ and e^g-j arc fixed the vectors e"^-^ and 6^2) are defined up to an arbitrary rotation in their own 
subspace. Transformations (4.46), (4.47) with _B = yield 



e 



— cos 8 e"-^^ + sin 6 6^2) , (4.63) 
e(^) = - sin 6 e"(i) + cos 9 e'^^) , (4.64) 

where Q is the rotation angle in the e"2)| plane. Vectors e"j^^ and 6^2) play a significant role in the discussion 

of polarized radiation. 



2. Propagation Laws for Reference Tetrad 

Discussion of the rotation of the polarization plane and the change of the Stokes parameters of electromagnetic 
radiation is not conceivable without an unambiguous definition of a local reference frame (tetrad) constructed along 
the light-ray geodesic. We start construction of this tetrad at the point of observation and prolongate it backward 
along the light-ray geodesic. By definition, the tetrad frames e"^^ and (/", n", to", m") are parallel transported along 
the light ray. The propagation equations for these vectors are thus obtained by applying the operator Dx = l'^'^ a of 
the parallel transport. For example, 

-i^^n.^'^\,)=^^ (4-65) 

where A is an affine parameter along the light ray. Using the definition of the Christoffel symbols (4.14) and changing 
over to the variable r with dx" /dr ~ k" + 0{h), one can recast equation (4.65) in the approximate form 



d_ 
d^ 



= \r {d.h^p - d^Kp) k^e^^ . (4.66) 



The law of propagation of the null vectors to" and m°' are 

^^+r^^//3TO^ = 0, (4.67) 

^+r^,?'^m'' = 0, (4.68) 

and the same laws are valid for n" and (sec equation (4.13)). Equations (4.66)-(4.68) are the main equations for 
the discussion of the rotation of the plane of polarization and variation of the Stokes parameters. 

3. Relativistic Description of Polarized Electromagnetic Radiation 

We consider propagation of plane electromagnetic waves. Such wave has an electromagnetic tensor Fap = J-ai3+0{e) 
defined in the first approximation by equation [102, 104] 

^ai3 = exp , (4.69) 

aap = * l[ami3] + * (4.70) 
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where $ is a complex scalar amplitude of the electromagnetic wave with two real components which are independent 
of each other in the most general case of incoherent radiation. In the proper frame of an observer with 4-velocity 
the components of the electric and magnetic field vectors are defined respectively as E" = —F^^up and H" = 
{-l/2y/^)e°''^"'^F^sui3. The electric field is a product of slowly-changing amplitude £a = —OLapu^ and fast-oscillating 
phase exponent 



Ea = £a CXp 



(4.71) 



The polarization properties of electromagnetic radiation consisting of an ensemble of waves with equal frequencies but 
different phases are defined in terms of the electric field measured by an observer. In the rest frame of an observer 
with 4-velocity , the intensity and polarization properties of the radiation are described in terms of the 'polarization 
tensor [1] 



Jafi ~< EaEfj > — < £aS/3 > , 



(4.72) 



where the angular brackets represent an average with respect to an ensemble of electromagnetic waves with randomly 
distributed and fast-oscillating phases. This averaging eliminates all fast-oscillating terms from Jaf3- One has to 
notice [1] that the polarization tensor Jaf3 is not fully symmetric since, in general, 8^ 7^ fa- Furthermore, the tensor 
Jafs satisfies the equalities Ja/su^ = Japl^ = which follows directly from its definition (4.72) and antisymmetry of 

The amplitude £a of the electric field can be decomposed in two components in the plane of polarization. Two tetrad 
vectors, m",m", form the circular polarization basis in this plane, and vectors e"]^^,e"2) form a linear polarization 
basis. The decomposition reads 



where 



1 ^ c 1 t: 



2 ' "2 
are left and right circularly-polarized components of the electric field, and 



(4.73) 
(4.74) 



(4.75) 



^(2) = ^ ($ - <&) 



(4.76) 



are linearly polarized components of the electric field, and uj — —lau" is the frequency of the electromagnetic wave. 

The set S = (/, Q, U, V) of electromagnetic Stokes parameters is defined with respect to two of the four vectors of 
the observer's tetrad e"^^ introduced in (4.48). We have [1, 105] 



I = Jaf3 

Q = JaP 

U = Jaf3 

V = iJap 



(1) (1) ^ C^r (2) 

„/3 



(l)"^ (1) 



(2)f= (2) 



p" pP -I- p" pf^ 

^ {ir (2) ^ ^ {"i-r (1) 



p/3 
• (1)^ (2) 



(2)6 (1) 



(4.77) 

(4.78) 
(4.79) 
(4.80) 



where / is intensity, Q and U - linear polarizations, and V - circular polarization of the electromagnetic wave. It is 
important to emphasize that though the Stokes parameters have four components, they do not form a 4-dimensional 
vector. This is because they do not behave like a vector under the Lorentz transformations [1]. 

Using equation (4.72), the Stokes parameters can be expressed in the standard way [1] in a linear polarization basis 



(l)'^(2) 



/ = <|£a)p + |£(2)P>, 

Q = <|£(1)P-|£(2)|2>, 

U = < £(i)£(2) + £'(1)£(2) > , 

V = i < £(i)£(2) - f (i)f(2) > 



(4.81) 
(4.82) 
(4.83) 
(4.84) 
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where = £"^6"^^ for n = 1, 2. Under the gauge subgroup of the httle group of l", the Stokes parameters remain 
invariant. However, for a constant rotation of angle Q in the (e"i),e"2)) polarization plane, one has 

I' = I , (4.85) 

Q' = Qcos2e + C/sin2e , (4.86) 

U' = C/cos2e-Qsin2e , (4.87) 

V = V . (4.88) 

This is what would be expected for a spin-1 field. That is, under a duality rotation of O = 7r/4, one linear polarization 
state turns into the other, while the circular polarization state remains the same. The transformation properties 
(4.85)-(4.88) of the Stokes parameters point out that the Stokes parameters Q, U represent the linearly polarized 
components, and V represents the circularly polarized component. 

The polarization vector P = (Pi, P2, -Ps) and the degree of polarization P = |P| of the electromagnetic radiation can 
be defined in terms of the normalized Stokes parameters by P = {Q/I, U/I, V/I). Any partially polarized wave may 
be thought of as an incoherent superposition of a completely polarized wave with the degree of polarization P and the 
polarization vector P, and a completely unpolarized wave with the degree of polarization 1 — P and zero polarization 
vector, P = 0, so that for arbitrary polarized radiation one has: (/, Q, U, V) = {PI, Q, U, V) + {I — PI, 0, 0, 0). For 
completely polarized waves, vector P describes the surface of the unit sphere introduced by Poincare [1]. The center 
of Poincare sphere corresponds to unpolarized radiation and the interior to partially polarized radiation. Orthogonally 
polarized waves represent any two conjugate points on the Poincare sphere. In particular, (Pi = ±1, P2 = 0, P3 = 0), 
and (Pi = 0, P2 = 0, P3 ~ ±1) represent orthogonally polarized waves corresponding to the linear and circular 
polarization bases, respectively. 



4- Propagation Laws for the Stokes Parameters 

Taking definition (4.70) of the electromagnetic tensor and accounting for the parallel transport of the null vectors 
Z", to", m" along the light ray and the laws of propagation of the electromagnetic tensor given by equations (4.15), 
yield the law of propagation of the complex scalar functions $ and $ 

— +79$ = 0, (4.89) 
aX 

^ + ?9$ = , (4.90) 

where A is an affine parameter along the ray and 'd = (1/2)Vq.P as given by equation (4.16). One notices that if A 
is a two-dimensional area of the cross section of the congruence of light rays lying in a surface of constant phase (j), 
then [2, 59, 102] 

-i^-2M = 0. (4.91) 
dX 

Thus, the product, .4|<I>p, remains constant along the congruence of light rays: 

Da iA\<^>\^) = . (4.92) 

This law of propagation for the product yl|(f>p corresponds to the conservation of photon's flux [2, 59]. 

The law of conservation of the number of photons propagating along the light ray, corresponds to the propagation 
law of vector ]$]^P. Indeed, taking covariant divergence of this quantity and making use of the equations (4.89), 
(4.90) along with definition (4.16) for the expansion of the bundle of light rays, yields [2, 59] 

Va d^l^Z") =0 . (4.93) 

This equation assumes that the scalar amplitude <f> of the electromagnetic wave must be interpreted in terms of the 
number density of photons in phase space, J\f, and the energy of one photon, i^photon = So;, as follows [2] 

/ A/- \ 1/2 

|$| ^ , (4.94) 

\ ^photon / 
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where the Planck reduced constant h = h/27r has been introduced to achieve consistency between the classical and 
quantum definitions of the energy of an electromagnetic wave. 

Each of the Stokes parameter is proportional to the square of frequency of light, ui = —u^la, as directly follows 
from equations (4.81)-(4.84) and (4.76). Therefore, the variation of the Stokes parameters along the light ray can be 
obtained directly from their definitions (4.77)-(4.80) along the light ray and making use of the laws of propagation 
(4.89)-(4.91). However, the set of the Stokes parameters S = (/, Q, U, V) is not directly observed in astronomy. Instead 
the set of four other polarization parameters {F^, Pi, P2, P3) is practically observed [59, 106]. Here P = {Pi, P2, P3) 
is the polarization vector as defined at the end of the preceding section, and F^^ is the specific flux of radiation (also 
known as the monochromatic flux of a light source [106]) entering a telescope from a given source. The specific 
flux is defined as an integral of the specific intensity (also known as the surface brightness [106]) of the radiation, 

= Ii^{lo,1), over the total solid angle (assumed <C 47r) subtended by the source on the observer's sky: 

I^{uj,l)dn{i) , (4.95) 



where I ~ (sin cos 0, sin sin 0, cos 0) is the unit vector in the direction of the radiation flow and dQ{l) = sint 
is the element of the solid angle formed by light rays from the source and measured in the observer's local Lorcntz 
frame. 

The specific intensity I^^ of radiation at a given frequency lj = 2ttv, flowing in a given direction, I, as measured in 
a specific local Lorcntz frame, is defined by 

^ ^(^"^^gy) . (4.96) 

(i(time)d(area)(i(frequency)d(solid angle) 

A simple calculation (see, for instance, the problem 5.10 in [107]), reveals that 

where h is the Planck's constant. The number density J\f is invariant along the light ray and does not change under 
the Lorentz transformation. Invariance of A/" is a consequence of the kinetic equation for photons (radiative transfer 
equation) which in the case of gravitational field and without any other scattering processes, assumes the following 
form [2] 

Equations (4.97), (4.98) tell us that the ratio /c^/cj^ is invariant along the light-ray trajectory, that is 

\ 3 

_luj_ _ I ^ 



(4.99) 



where wq, ^ are frequency of light at the point of emission and observation respectively, Iu}o is the surface brightness 
of the source of light at the point of emission, and /^^ is the surface brightness of the source of light at the point of 
observation. 

Equations (4.89)-(4.99) make it evident that in the geometric optics approximation the gravitational field does not 
mix up the linear and circular polarizations of the electromagnetic radiation but can change its surface brightness 
due to gravitational (and Dopplcr) shift of the light frequency caused by the time-dependent part of the gravitational 
field of the isolated system emitting gravitational waves. Furthermore, the monochromatic flux from the source of 
radiation changes due to the distortion of the domain of integration in equation (4.95) caused by the gravitational 
light-bending effect. Taking into account that the gravitationally-unperturbed solid angle dVL{k) = smOdOdcf), and 
introducing the Jacobian, J{9, <j)), of transformation between the spherical coordinates {0,(j)) and {9,(j)) at the point of 
observation, one obtains that the measured monochromatic fiux is 

F^^ J d^ J I^,{e,(j))Ji9,q^)iiu/u;ofsm9d9 . (4.100) 

Equation (4.100) tells us that the monochromatic flux of the source of light can vary due to: 

1. the gravitational Dopplcr shift of the electromagnetic frequency of light when it travels from the point of emission 
to the point of observation; 
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2. the change in the solid angle at the point of observation caused by the gravitational light deflection. The 
"magnification" matrix is the Jacobian of the transformation of the null directions on the celestial sphere 
generated by the bending of the light-ray trajectories by the gravitational field of the isolated system. 

Orientation of two components, Pi and P2, of the polarization vector P at the point of observation differs from 
that taken at the point of emission of the electromagnetic signal due to the parallel transport of the reference tetrad, 
which changes its orientation with respect to the reference tetrad at infinity as it propagates along the light ray. 
This can be seen after taking into account equations (4.63), (4.64) where the rotational angle @ is determined as 
a solution of the parallel transport equation (4.66) for the reference tetrad. Thus, the gravitational field changes 
the tilt angle of the polarization ellipse as light propagates along the light-ray trajectory. This effect was predicted 
by Skrotskii [108, 109] and discussed independently in [110, 111]. Its observation will play a significant role in a 
foreseeable future for detection of gravitational waves of cosmological origin by CMBR-radiometry space missions 
[112]. The third component of the polarization vector, P3, remains the same along the light-ray trajectory because 
it represents the circularly polarized component of the radiation and is not affected by the rotation of the reference 
tetrad as it propagates along the light-ray path. 



5. THE MATHEMATICAL TECHNIQUE FOR ANALYTIC INTEGRATION OF THE LIGHT-RAY 

EQUATIONS 

This section provides mathematical technique for performing integration of equations of propagation of various 
characteristics of electromagnetic wave from the point of emission of light Xq to the point of its observation x. 
The basic element, which is to be integrated, is the metric tensor perturbation hap(t,x) generated by the localized 
astronomical system and taken at the point lying on the light-ray trajectory. The metric tensor depends on the 
multipole moments taken at the retarded instant of time s = t — r divided by the distance from the system r, and its 
time and spatial derivatives. These derivatives are handled by the method shown in this section. The main integrals 
are taken from the quantity F{s)/r, wherein F{s) denotes a multipole moment of the gravitating system depending 
on the retarded time s = t — r. Calculation of the integrals and their derivatives is slightly different for stationary 
and non-stationary components of the metric tensor and are treated in the next two sections. 



5.1. Light-ray Integrals from Monopole and Dipole 

The monopole and dipole parts of the metric tensor (3.5)^(3.8) are formed by terms being proportional to the 
mass A4 and spin 5* of the isolated system. These terms as they appear in the solution of the light-ray equations, 
are given by integrals [A4(t — 7')/r][~^l, [Ai{t — r)/r]t~^l and [S^(t — r)/r]^~^\ Wit — r)/r][~^l respectively. In this 
section we shall assume for simplicity that mass A4 and spin 5* of the isolated system arc constant during the time of 
propagation of light from the source of light to observer. The assumption about constancy of the mass Ai and spin 
iS* of the isolated system is valid as long as one neglects the energy emitted in the form of the gravitational waves 
by the isolated system under consideration. For light sources at the edge of our visible universe the characteristic 
time for emission of gravitational waves by an isolated system can be comparable with time the light takes to travel 
from the point of emission to observer. In this case the time evolution of the mass-monopole and spin-dipole is to be 
taken into account for correct calculation of the perturbations in propagation of light ray (see section 6.1 for more 
detail). Mass and spin of the isolated system can also change due to catastrophic disruption of the isolated system, 
for example, as a result of a supernova explosion. Specific details of how this process affects the light-ray propagation 
are not considered in the present paper but perhaps are worthwhile to studying. 

In the case when mass and spin are constant, the integrals that we need to carry out are reduced to [l/r]["^l and 
[l/r][~^l which are formally divergent at the lower limit of integration at the past null infinity when time t — > —00. 
However, one must bear in mind that these integrals do not enter equations of light-ray geodesies (4.37) alone but 
appear in these equations after taking at least one partial derivative with respect to either or t parameters. This 
differentiation effectively eliminates divergent parts of the integrals from the final result. Hence, in what follows, we 
drop out the formally divergent terms so that the integrals under discussion assume the following form 
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dT 
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(5.1) 
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where T = t — t* and we used equation (4.22) explicitly while expressing the distance r = VNir) as a function of 
time T. The constant distance was introduced to make the argument of the logarithmic function dimensionless. 
This constant is not important for calculations as it always cancel out in final formulas. However, in the case of 
gravitational lensing it is convenient to identify the scale constant with the radius of the Einstein ring [59-61] 

where r, tq are distances from the isolated system (the deflector of light) to observer and the source of light, and 
R = T + tq. The radius of the Einstein ring is a characteristic distance separating naturally the case of weak 
gravitational lensing (d > Te) from the strong one {d < r^). The Einstein ring as visible in the sky, has an angular 
size that is given by 



Te / AGM To 



r \ c? Rr 



1/2 

(5.4) 



The angular Einstein radius 9^ defines the angular scale for a lensing situation. For cosmology, typical mass M ~ 
10^^ Mq and distances r, tq, and R are of order 1 gigaparsec. Consequently, the angular Einstein radius 0^ is of 
order one arcsecond and the linear radius r-^ is of order 1 kiloparsec. A star within our galaxy having Af ~ Mq and 
distances r, ro, and R of order 10 kiloparsec has an angular Einstein radius 6^ of order a milli-arcsecond and the 
linear radius of order 10 astronomical units (10 AU). 



5.2. Light-ray Integrals from Quadrupole and Higher-Order Multipoles 

Time-dependent terms in the metric tensor (3.5) ~ (3.8) result from the multipole moments which can be either 
periodic (a binary system) or aperiodic (a supernova explosion) functions of time. The most straightforward way to 
calculate the impact of the gravitational field of such a source on the propagation of light is to decompose its multipole 
moments F{s) {Xl{s),Sl{s)) in the Fourier series [83, 113] 

F{t -r)^ —= / F((I;)e*^(*-'')(i(i , (5.5) 

V27r J -no 



substitute this decomposition to the light-ray propagation equations and integrate term by term. This method makes 
an impression that in order to obtain the final result of the integration the (complex-valued) Fourier image FiCo) of 
the multipole moments of the isolated system must be specified explicitly, otherwise the convolution of the integrated 
Fourier series will be impossible. However, this is not true, at least in general relativity, and, as we shall show later, 
the explicit structure of F{u)) is irrelevant for subsequent general-relativistic calculations. This is explained, if one 
recollects that in general relativity gravitational field propagates with the same speed as light. In alternative theories 
of gravity the speed of gravity and light may be different, so it can lead to appearance of terms proportional to the 
difference between the speed of gravity and light that will drastically complicates calculations which will require the 
Fourier images of gravitational-wave sources as explicit functions of w [114]. All such terms are, however, cancelled 
out in general relativity exactly, thus making calculations manageable and applicable for any source of gravitational 
waves. 

We notice that the integration of the light-ray equations is effectively reduced to the calculation of only two types 
of integrals along the light ray: [F{t — r)/r][~^l and [F{t — 7')/r][~^l, where F{t — r) denotes any type of the time- 
dependent multipole moments of the gravitational field of the localized astronomical system. These integrals can be 
performed after introducing a new variable [83] 



y 



t* = T- r(r) = r - ^d^ + , (5.6) 



which has a unique physical meaning of the retarded interval of time between two space-time events: position of 
photon = {t,x) on the light ray and the center of mass of the isolated system z" = (y, 0). It is important to 
realize that the equation (5.6) is a retarded solution of the null cone equation in flat space-time 

r;„;3(z"-z")(x'3-z^')=0, (5.7) 



describing propagation of gravity from the isolated system to the photon. This retardation of gravity effect has its 
origin in the time argument of the solution (3.5)-(3.8) of the linearized Einstein equations and is due to the finite 
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speed of propagation of gravity. The retardation of gravity effect was predicted in [115] and experimentally tested 
with the precision 20% in the jovian light-deflection experiment in September 2002 [116]. This replacement of the 
time argument from r to the retarded time y = r — r(r) allows us to perform integration of equations of light-ray 
geodesies completely without making specific assumptions about time dependence of the multipole moments. It is 
worth noting that while the parameter r runs from — c» to -l-oo, the retarded time y runs from — oo to 0; that is, y is 
always negative, j/ < 0. 

Equation (5.6) leads to other useful transformations 



y 



IcP + y^ 



(5.8) 



22/ ' ^ 2 y ' 

Making use of the new variable y and relationships (5.8) the integrals under discussion can be displayed as follows 
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(5.10) 



— oo — oo 



-oo — oo 



where C, rj are new variables of the integration replacing integration along the light-ray trajectory by that along a 
null characteristic of the gravitational field, and t* is the time of the closest approach of photon to the origin of the 
coordinate system. The time t* has no special physical meaning in general case and appears in calculations as an 
auxiliary (constant) parameter which vanishes in the final result. It can be used as an approximation of the retarded 
time s ~ t — r in the case of a small impact parameter d of the light ray to the isolated astronomical system, that 
is s ~ t* -f- 0{d^ /cr). Notice that for observer at spatial infinity r = oo and the retarded time s ~ t* which justifies 
dependence of the total angle of the deflection of light on. otherwise unphysical, value of t* adopted, for example, in 
[82]. 

A remarkable property of the integrals in the right side of equations (5.9), (5.10) is that they depend on the 
parameters ^ and t only through either the upper limit of the integrals which is the variable y = t — \^(P + t^, or 
the square of the impact parameter, — ■ ^) standing in front of the second integral in the right side of equation 
(5.10). For this reason, differentiation of the integrals in the left part of equations (5.9), (5.10) with respect to either 
^* or T will effectively eliminate the integration along the light ray trajectory. For example. 



F{t-r) 
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(5.11) 



Similar calculation can be easily performed in case of differentiation of integral [F{t — r)/r][ ^1 with respect to t. It 
results in 



dr 



F{t - r) 
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(5.12) 



as it could be expected because the partial differentiation with respect to r keeps ^* and t* fixed and, hence, is 
equivalent to taking a total time derivative with respect to t along the light ray. 

Since all terms in the solution of the light geodesic equation are represented as derivatives from the integrals 
[F{t — r)/r][~^l and [F{t — r)/r]I~^l with respect to the parameters C and/or r, it is clear from equations (5.11), 
(5.12) that the solution will not contain any single integral [F{t — r)/r][~^l at all; only integrands of this integral will 
appear. We notice that taking first and second derivatives of the double integral \F{t — r)/r][~^l do not eliminate the 
integration along the light ray trajectory 
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However, taking one more (a third) derivative eliminates all integrals from equation (5.14) completely. More specifi- 
cally, we have 
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and making use of equations (5.11), (5.14) for explicit calculation of partial derivatives in the right part of equation 
(5.15) proves that all integrations disappear, that is 
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where F(t — r) = dtF{t — r). The same kind of reasoning works for the third derivatives with respect to the parameter 
T and to the mixed derivatives taken with respect to both ^* and r. 

We shall obtain solution of the light geodesic in terms of STF derivatives with respect to the impact parameter 
^* of the light ray acting on the single and double integrals having a symbolic form d^ai...ak>[F{t ~ 7')/r][~^l and 
d<ai...ak>[F{t ~ ?')/r][~^l. Explicit expression for this STF derivative can be obtained directly by applying differ- 
entiation rules shown in equations (5.13) - (5.15). In what follows, we shall see that solutions of equations of light 
rays will always contain three or more derivatives acting on the integrals from the higher-order multipole moments 
(multipolarity I > 2) which have the same structure as that shown in equations (5.9), (5.10). It means that the 
final result of the integration of the light-ray propagation equations depending on the higher-order multipoles of the 
localized astronomical system can be expressed completely in terms of these multipoles taken at the retarded instant 
of time s = t — r. In other words, motion of photon in general relativity does not depend on the past history of 
its propagation and the effects caused by temporal variations of the gravitational field (gravitational waves) are not 
accumulating [117]. On the other hand, the past history of the isolated system can affect the propagation of the light 
ray through the tails of the gravitational waves contributing to the multipole moments of the system as given by 
equations (3.12), (3.13). 



6. RELATIVISTIC PERTURBATIONS OF THE LIGHT- RAY PROPAGATION 



6.1. Relativistic Perturbation of the Electromagnetic- Wave Eikonal 

Perturbation of eikonal (phase) of electromagnetic signal propagating from the point Xq to the point x are obtained 
by solving equation (4.35). This solution is found by integrating the metric tensor perturbation along the light-ray 
trajectory and can be written down in the linearized approximation as an algebraic sum of three separate terms 
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where ip^cf represents the gauge-dependent part of the eikonal, and iI'(m), V'(s) are eikonal's perturbations caused by 
the mass and spin multipole moments correspondingly. Their exphcit expressions are as fohows 
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where H{p — q) is the Heaviside function defined by the expression (f .18), C{p, q) is the polynomial coefRcient given 
by equation (1.17) and F'-P^{t - r) denotes dPF{t - r)/dtP where r is considered as constant. 

We note that the gauge-dependent part of the eikonal, '0(g), contains combination of terms fc'iys' — ip^ defined by 
equations (C.7)-(C.9) which can be, in principle, eliminated with an appropriate choice of the gravitational field gauge 
functions uP and w*. However, such a procedure will introduce a reference frame in the sky with a coordinate grid 
being sensitive to the direction to a source of light rays; that is, to the unit vector fc*. The coordinate frame obtained 
in this way will have the direction-dependent distortions leading to its trifling practical usage. For this reason we do 
not recommend the elimination of functions and tp^ from equation (6.2) and give preference to the ADM-harmonic 
coordinate system which admits a much simpler and unique treatment of observable relativistic effects. Thus, we 
leave the functions ip^ and in equation (6.2) where the gauge functions vlP and w"^ are defined by formulas (C.2), 
(C.3). 

Expressions (6.3), (6.4) for the mass- and spin-multipoles of the eikonal contains integrals. However, one can easily 
prove that the integrals from all high-order multipoles are eliminated. Indeed, scrutiny inspection of equations (6.3), 
(6.4) elucidates that all the integrals from the high-order multipole moments enter the equation in combination with 
at least one derivative with respect to the impact parameter vector Hence, the differentiation rule (5.11) is applied 
which eliminates the integrals from the high-order multipoles in the eikonal. 

The integral, which remains unperformed, is that from the mass monopole in equation (6.3). It can be taken by 
parts as follows 
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If one assumes that the mass M. is constant, then the second term in the right side of equation (6.5) vanishes and 
the eikonal does not contain any integral dependence on the past history of the light propagation. This assumption 
is usually implied, for example, in the theory of gravitational lensing [59-61] and other numerous applications of the 
relativistic theory of light propagation. Here, we extend our approach to take into account the case of 7^ 0. In the 
most general case the lost of energy by the isolated system caused by emission of gravitational waves is given by [1, 2] 
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where Xj^J' represents a third time derivative from the mass quadrupole moment, and terms of order O [c 
contribution of higher-order mass and spin multipoles. Accounting for the contribution of the gravitational waves 
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given by equation (6.6) yields 
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It is instructive to evaluate contribution of the second term in the right side of equation (6.5) in the case of hght 
propagating in gravitational field of a binary system consisting of two stars with masses mi and m2 orbiting each 
other along a circular orbit. The total mass M of the system is defined in accordance with equation (3.10) which 
takes into account the first post-Newtonian correction 



M{t) = mi + m2 
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where ^ = mim2/A^ is the reduced mass of the system and a = a(t) is the orbital radius of the system. Assuming 
that the lost of the orbital energy is due to the emission of gravitational waves from the system in accordance with 
the quadrupolc formula approximation [1, 2], the time evolution of the orbital radius is defined by equation [1, 2] 
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It has a simple solution [2] 
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where t = t — t* , t* is the time of the closest approach of light ray to the binary system, gq = a{t*) is the orbital 
radius of the binary system at the time of the closest approach, and T is the spiral time of the binary. The lost of 
the total mass due to the emission of the gravitational waves is 
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Substituting equation (6.12) to the right side of equation (6.5) and performing integration yields 



2 c^a{s) 



ajs) - ap 
a(s) + ao 



2 arctan 



a(s) 
ao 



where a(s) = a{t — r) denotes the orbital radius of the binary system taken at the retarded time s = t 
all terms in equations (6.5), (6.8), (6.13) together yields 
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Putting 



(6.14) 



where the first term in the right side of this equation is the standard Shapiro time delay in the gravitational field of 
the binary system with constant total mass mi + m2, and the second term represents relativistic correction due to 
the emission of gravitational waves by the system causing the overall loss of its orbital energy. 

Eikonal describes propagation of a wave front of the electromagnetic wave. Light rays are orthogonal to the wave 
front and their trajectories can be easily calculated as soon as the eikonal is known. In the present paper we do 
not use this technique and obtain solution for the light rays directly from the light-ray geodesic equations that gives 
identical results. 



6.2. Relativistic Perturbation of the Light-Particle Velocity 

Integration of the light-ray propagation equations (4.40) - (4.43) is fairly straightforward. Performing one integra- 
tion of these equations with respect to time yields 

i'(r) = fc' + S'(r,0 , (6.15) 

5''(r, ^) = S» (r, + S» (r, |) + (r, |) , (6.16) 
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where the relativistic perturbations of photon's trajectory are given by 

(r, = dr [{if' - fcV°) + iw' ~ fc'w°)] , (6.17) 
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Here H{p— q) is a Heaviside function defined by the expression (1.18) and Ci are the polynomial coefficients (1.17). 
The gauge functions are given in appendix C. Mass monopole and spin dipole terms are written down in equations 
(6.18), (6.19) in symbolic form and after taking derivatives are simplified 
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Remaining integrals shown in the right side of equations (6.18), (6.19) are convenient for presentation of the result of 
integration in the symbolic form. The integrals are actually eliminated because they are always appear in combination 
with, at least, one operator of derivative di with respect to the unperturbed value of the impact parameter of the 
light ray. The derivative operator acts on the integrals in accordance with equation (5.11) converting the integrals 
into fimctions of the retarded time y = t — rir) 
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We conclude that at each point of the wave front of the electromagnetic wave the relativistic perturbation of the 
direction of propagation of light ray (wave vector) caused by the time-dependent gravitational field of the isolated 
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system depends only on the value of its multipole moments taken at the retarded instant of time and it does not 
depend on the past history of the light propagation. 

6.3. Perturbation of the Light-Ray Trajectory 

Integration of equation (6.15) with respect to time yields the relativistic perturbation of the trajectory of the light 
ray 
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Herein the term 
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is the gauge-dependent part of the trajectory's perturbation, and the physically meaningful perturbations due to the 
mass and spin multipoles 
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Here i/(p — q) is again a Heaviside function defined by equation (1.18) and Ci are the polynomial coefficients (1.17). 
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Relativistic perturbations (6.27), (6.28) contain two types of integrals along the light ray - [F{t — r)/r][~^l and 
[F{t — r)/r][~^l. In section 5 we have shown that taking a derivative from the integral [F{t — r)/r][~^l effectively 
reduces the integral to the ordinary function of the retarded time s = t — r with no dependence on the past history 
of the propagation of light ray. In order to eliminate second order integrals [F{t — r)/r]^~^^ one needs to take three 
or more derivatives. One can notice that in (6.27), (6.28) there are some terms in which the number of derivatives 
is not sufficient to eliminate the integrals. Nevertheless all such terms in (6.27), (6.28) are cancelled out because 
the numerical factors standing in front of these terms get nullified. As an example let us examine equation (6.27) 



for (t, ^). In this expression the integral '^\t — r)/r 



(Af) 



is differentiated I — p + 1 times with respect to 



the impact parameter (derivative da) and q times with respect to the parameter r (derivative dr)- Expression 

9<ai...a,_p+i>3^ '"{t — r)/r can contain the integral dependence if either (a) q = 0, p e (/ — 1, or (b) 9 = 1, 

p ^ I. In any case either factor I — (p — q)/l or that 1 — {p — q)/ {I — 1) will be zero, thus, annihilating the integrals 
from the mass multipoles in equation (6.27). Similar consideration can be done for the spin-dependent perturbation 
(r, ^) which shows that there is no integral dependence of the spin multipolc moments on the past history of the 

(S) 

light propagation. 

We emphasize, however, that the integral dependence of the light-ray perturbation on the past history remains in 
equation (6.26) which contains the gauge functions ip" and w" used later in this paper for interpretation of observable 
effects of the gravitational waves. The past-history dependence of the light-ray trajectory comes into play also 

[—2] [—21 

through the integrals from mass monopole and spin-dipole terms, [A4/r] and [iS/r] , if mass and/or spin are not 
conserved and change as time passes. The non-conservation can be caused, for example, by emission of gravitational 
waves carrying away the orbital energy and angular momentum of the system. The past-history contribution of these 
integrals can be calculated similarly to eikonal perturbation (see equations (6.6)-(6.14)) but we do not carry out this 
calculation in the present paper. 

Solution of the light-ray equation with the initial-boundary conditions (4.17) depends on the unit vector k defining 
direction of the light-ray propagation extrapolated backward in time to the past null infinity. In real practice the 
light ray is emitted at the point Xq where the source of light is located at time to, and it arrives to observer at time 
t to the point x separated from the source of light by finite coordinate distance R = |a; — a;o|. Therefore, solution of 
the light-ray equations must be expressed in this case in terms of the integrals of the boundary-value problem. It is 
formulated in terms of initial Xq and final x positions of the photon 

x{t) ~ X , x{to) = Xq , (6.29) 

and assumes that the solution of the light-ray equation depends on the unit Euclidean vector 

A- = (6.30) 

\X - Xq\ 

that defines the direction from the observer towards the source of light and may be interpreted as a coordinate 
direction in the Euclidean space. 

In what follows it is convenient to make use of the astronomical coordinates a; = x' = {x^, a;^, x^) based on a triad 
of the unit vectors (/q, Jq, Kq) (see Figure 2). Vector Kq points from the observer towards the deflector, and vectors 
Iq and Jo lie in the plane of the sky defined as the plane being orthogonal to vector Ko- The unit vector Jq is directed 
to the east, and that Jq points towards the north celestial pole. The origin of the coordinate system is chosen to lie 
at the barycenter of the isolated system which emits gravitational waves. 

Another reference frame based on a triad of the unit vectors {I, J, K) rotated with respect to vectors (/q, Jq, Kq) 
is useful as well. Vector K points from the observer towards the source of light, and vectors / and J lie in the plane 
of the sky defined to be orthogonal to vector K . We emphasize that this plane is different from the plane of the sky 
being orthogonal to vector Kq- This is because the plane of the sky is actually a tangent plane to the unit sphere, 
and vectors K and Ko point in different directions on it. Mutual orientation of one triad with respect to another one 
is determined by the following orthogonal transformation 

Jo = J cos -I- J sin 51 , (6.31) 
Jo = —Icos9smil. + Jcos9cosft + Ksm9, (6.32) 
Ko = / sin 6' sin fl — J sin 6* cos + JtT cos , (6.33) 

where the rotational angles CI and 9 are constant. 

It is rather straightforward to obtain solution of the boundary value problem for light propagation in terms of the 
unit vector K. All what one needs is to convert the unit vector k to K written in terms of spatial coordinates of the 
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points of emission. Xq, and observation, of the light ray. From formula (6.24) one has 
where the relativistic correction (r, ^) to the vector A"' is defined as follows 
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and P\ 



^ „^ ,„ is the operator of projection on the plane being orthogonal to vector fc*. Equation (6.35) 
contains in its denominator the distance R — \x ~ Xq\ between observer and source of light which can be very large 
making impression that the relativistic correction /3'(t, ^) is negligibly small. However, the difference AS^(t, ^) in 
the numerator of (6.35) is proportional either to the distance r = |a;| between observer and the isolated system or 
to the distance rp = |a;o| between the source of light and the isolated system. Both distances can be comparable 
with R so that the relativistic correction /3*(t, ^) can be observable and must be taken into account, in general, for 
calculation of light deflection in the cases of finite distances of observer and/or source of light from the localized 
source of gravitational waves. Only in the case where observer and source of light reside at extremely large distances 
on opposite sides of the source of gravitational waves can the relativistic correction /?' be neglected. 



7. OBSERVABLE RELATIVISTIC EFFECTS OF THE MULTIPOLAR GRAVITATIONAL FIELDS 



In this section we derive general expressions for four observable relativistic effects - time delay, light bending, 
frequency shift, and rotation of the polarization plane. 



7.1. Gravitational Time Delay of Light 



Relativistic time delay can be obtained cither directly from the expression (6.24) or from the electromagnetic 
eikonal (4.32), (6.1) by observing that the eikonal is constant not only on the null hypersurface of the wave front of 
the electromagnetic wave but also along the light rays. Both derivations lead to the same result 

t-<o=|a;-a;o| + A(r,To) , (7.1) 
A(r, ro) - A (t, tq) + A (t, tq) + A (t, tq) , (7.2) 
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where Xq = (to.xo) are coordinates of the point of emission of light and x" = {t,x) are coordinates of the point of 
observation, A (r, tq). A (r, tq) and A (r, tq) are functions describing the delay of the electromagnetic signal due to 
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the gauge-freedom, mass and spin multipoles of the gravitational field of the isolated system correspondingly. These 
functions are expressed as follows 
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(7.3) 
(7.4) 
(7.5) 



We also notice a relationship between the relativistic perturbations of the eikonal and the light-ray trajectory 

^(r,0 = -fc,S'(r,0, (7.6) 

where the eikonal perturbation ij) reads off the equations (6.1)-(6.4). From equations (7.3)~(7.5) one can infer that in 
the linear approximation the functions describing the time delay are just the projections of the coordinate perturbation 
of the light-ray trajectory onto the unperturbed direction A:* from the source of light to observer. 

It is worth to remark that equation (7.1) defines the time delay effect with respect to the global ADM-harmonic 
reference frame. In order to convert it to observable proper time, we assume for simplicity that the observer is in a 
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state of free fall and moves with veloeity with respeet to the reference frame of the isolated system. Transformation 
from the ADM-harmonic coordinate time t to the proper time T is made with the help of the standard formula [1, 2] 

dT^ = -gafsdx^dx^ . (7.7) 

Substituting the metric tensor expansion ga/s = j/q/j + /iq/3, where hap is given by equations (3.21)-(3.24), to equation 

(7.7) yields 

^ = ^1 - - hoo{l + V2) - 2ho^V^ - h^V^Vi . (7.8) 

In the most simple case, when observer is at rest {V = 0) with respect to the ADM-harmonic reference frame equation 

(7.8) is drastically simplified and depends only on /iqo component of the metric tensor. Implementation of formula 
(3.21) for hoo and subsequent integration of (7.8) with respect to time yields then 



T = (^1-— j (t-U) , (7.9) 

where ti is the initial epoch of observation. Another simple case of equation (7.8) is obtained for observer located at 
the distance r so large that one can neglect /iqo f-nd ft-oi quasi-static perturbations of the metric tensor. Then, the 
difference between the observer's proper time T and coordinate time t is 

^ = ^l-V^-hj;^V^V^ , (7.10) 

leading in the case of a small velocity to 

T = Vl-V^{t-U)-V'V^ hf^dt. (7.11) 

The last term in the right side of this equation describes the effect of the gravitational waves on the rate of the proper 
time of observer which may be important under some specific circumstances but is enormously small in real practice 
and; as a rule, can be neglected. 



7.2. Gravitational Deflection of Light 

The coordinate direction to the source of light measured at the point of observation x is defined by the four- vector 
of photon P. This vector depends on the energy (frequency) of photon which is irrelevant for the present section. 
Normalization of I" to its frequency brings about four- vector = (l,p*) where = — x', which is proportional to 
P but is only direction-dependent. Spatial part of this vector can be decomposed around direction k' of unperturbed 
photon's trajectory 

p' = -k' - , (7.12) 

where overdot denotes derivative with respect to coordinate time t and the minus sign indicates that the tangent 
vector p^ directs from observer to the source of light. 

The coordinate direction is not a directly observable quantity as it is defined with respect to coordinate grid 
on the curved space-time manifold. A real observable vector towards the source of light, s" = (l,s*), is defined 
with respect to the local inertial frame co-moving with the observer. In this frame s* = —dX^/dT, where T is the 
observer's proper time and X' are spatial coordinates of the local inertial frame with observer being located at its 
origin. We shall assume for simplicity that observer is at rest with respect to the global ADM-harmonic coordinate 
system {t,x^). The case of observer moving with respect to the ADM-harmonic system with velocity may be 
considered after completing the additional Lorentz transformation which is pretty straightforward procedure so that 
we do not discuss it from now on. In case of the static observer the infinitesimal transformation from the global 
ADM-harmonic coordinates (t,.i;*) to the local coordinates (T, AT*) is given by the formulas 

dT = A°, dt + A° dx^ , dX' = A*o dt + A^ dx' , (7.13) 

where the matrix of transformation A^ is defined by the requirements of orthonormality 



(7.14) 
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In particular, the orthonormality condition (7.14) assumes that spatial angles and lengths at the point of observations 
are measured with the Euclidean metric S^j. Because the vector s" is isotropic (null) with respect to the Minkowski 
metric rjafj, we conclude that the Euclidean length |s| of the vector s* is equal to 1. Indeed, one has 

7?«/3s"s'^ = "l + s2 = 0. (7.15) 

Hence, |s| = 1. 

In the linear approximation with respect to the universal gravitational constant G, the matrix of the transformation 
is as follows [118] 



A°o = 1 - ^hon{t,x) , 



A", = -ha,{t,x) , 
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(7.16) 



Using transformation (7.13) we obtain relationship between the observable vector s* and the coordinate direction p' 
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In the linear approximation this takes the form 
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Remembering that vector \s\ — 1, wc find the Euclidean norm of the vector from the relationship 



1 



IpI = 1 - /loo + hojP^ - ^hj^py 



which brings equation (7.18) to the form 



(7.17) 



(7.18) 



(7.19) 



(7.20) 



where P'^ = 5"^^ — k^k^ , and the Euclidean unit vector m* = pVIpI- 

Let us now denote by a* the dimensionless vector describing the total angle of deflection of the light ray measured 
at the point of observation, and calculated with respect to vector fc' given at past null infinity. It is defined according 
to the relationship [119] 



/{t,0 - fcMfc-H(r,0]-S*(r,0, 



or 



(7.21) 
(7.22) 
(7.23) 



a'(r,0 - -P^S^(r,0. 
As a consequence of the definitions (7.12) and (7.22) we conclude that 

TO* = -k' + a\T,^) . 

Accounting for expressions (7.20), (7.23), and (6.34) we obtain for the observed direction to the source of light 

(r, ^) = + a' (r, ^ ) + (r, ^ ) + f (r, ^ ) , (7.24) 

where the unit vector is given by equation (6.30), relativistic correction /3* is defined by equation (6.35), and the 
perturbation 



1 



7Xr,0 = ~-P'^k'^hl^it,x) 



(7.25) 
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determines deformation of the focal coordinates of observer with respect to the gfobal ADM-harmonic frame caused 
by the transverse-tracekss part of the gravitational field of the isolated system at the point of observation. 

If two sources of light (quasars, stars, etc.) are observed along the directions Si and S2 corresponding to two different 
light rays passing near the isolated system with impact parameters, ^1 and ^2, the measured angle ^' between them 
in the focal inertial frame is 



cos 5* 



(7.26) 



where the dot between the two vectors denotes the usual Euclidean scalar product. It is worth emphasizing that 
the observed direction s* to the source of light includes relativistic deflection of the light ray in the form of three 
perturbations. Two of them, a' and 7*, depend only on quantities taken at the point of observation but according 
to equation (6.35), is also sensitive to the strength of the gravitational field taken at the point of emission of light. 
This remark reveals that according to relation (7.24) a single gravitational wave signal may cause different angu- 
lar displacements and/or time delays for different sources of light located at different distances from the source of 
gravitational waves even if the directions to the light sources are the same. 

Without going into further details of the observational procedure we give an explicit expression for the total angle 
of the light deflection a*. We have 
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7.3. Gravitational Shift of Frequency 



The exact calculation of the gravitational shift of electromagnetic frequency between emitted and observed photons 
plays a crucial role for the adequate interpretation of measurements of radial velocities of astronomical objects [120], 
anisotropy of electromagnetic cosmic microwave background radiation (CMBR), and other spectral astronomical 
investigations. In the last several years, for instance, radial velocity measuring technique has reached unprecedented 
accuracy and is approaching to the precision of about 10 cm/sec [121]. In the near future there is a hope to improve the 
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accuracy up to 1 cm/sec [122] when measurement of the post-Newtonian relativistic effects in spectroscopic binaries 
and/or multiple star systems will be possible [120]. 

Let a source of light move with respect to the ADM-harmonic coordinate frame (t, a;*) with velocity Vq = Xo{to) 
and emit electromagnetic radiation with frequency vq = l/((5To), where to and Tq are the coordinate and proper 
time of the source of light, respectively. We denote by = 1/{ST) the observed frequency of the electromagnetic 
radiation measured at the point of observation by an observer moving with velocity V ~ x{t) with respect to the 
ADM-harmonic coordinate frame. In the geometric optics approximation we can consider the increments of the proper 
time, STq and 6T, as infinitesimally small which allows us to operate with them as with differentials [123]. Time delay 
equation (7.1) can be considered as an implicit function of the emission time to = to{t) having the time of observation 
t as its argument. Because the coordinate and proper time at the points of emission of light and its observation are 
connected through the metric tensor we conclude that the observed gravitational shift of frequency 1 + 2 = iz/i/q can 
be defined through the consecutive differentiation of the proper time of the source of light, Tq, with respect to the 
proper time of the observer, T, [123, 124] 

dTo (ITq cHq dt 

Synge calls the relationship (7.31) as the Doppler effect in terms of frequency (see [124], page 122). It is fully 
consistent with definition of the Doppler shift in terms of energy (see [124], page 231) when one compares the energy 
of photon at the points of emission and observation of light. The Doppler shift in terms of energy is given by 

1 + z = — = g , (7.32) 

where Mq, u" and Iq, I" are the 4- velocities of the source of light and observer and the 4-momenta of photon at 
the points of emission and observation respectively. It is quite easy to see that both mentioned formulations of the 
Doppler shift effect are equivalent. Indeed, taking into account that = dx^/dT and la — d(p/dx°', where (p is the 
phase of the electromagnetic wave (eikonal), we obtain rjapu^l^ = dip/dT. Thus, equation (7.32) yields 

l + ^ = ;^§^. (7.33) 

dipo dl 

The phase of electromagnetic wave remains constant along the light ray trajectory. For this reason, dLp/dLpi^ = 1 
and, hence, equation (7.33) is reduced to equation (7.31) as expected on the ground of physical intuition. Detailed 
comparison of the two definitions of the Doppler shift and the proof of their identity in general theory of relativity is 
thoroughly discussed in our paper [125]. 

We emphasize that in equation (7.31) the time derivative 



dTo_ 
dto 



l-Vo^-{l + V„^)hooito,xo) - 2V^hH{to. xo) - V^V^hf^{to, xo)] , (7.34) 



is taken at the time tQ at point of emission of light Xq, and the time derivative 

^ = [l-V^-il + V')hoo{t,x)~2V'ho,it,x)~V^V^hf;^it,x)]-'^\ (7.35) 

is calculated at the time of observation t at the position of observer. 

The time derivative dto/dt along the light-ray trajectory is calculated from the time delay equation (7.1) where 
we have to take into account that function A(t, To) depends on times to and t indirectly through the retarded times 
So ~ to ~ ro and s ~ t — r that are arguments of the multipole moments of the isolated system and wherein 
rg = |a;o| — \x{to)\, r = |a;| = \x{t)\ are functions of time to and t respectively. Function A(t, tq) also depends on the 
time of the closest approach of light ray, t* , through variables t = t — t* , tq = to — t*, and on the unit vector k. Both 
t* and k should be considered as parameters depending on to and t because of the relative motion of observer with 
respect to the source of light which causes variation in the relative position of the source of light and observer and, 
consequently, to the corresponding change in the parameters characterizing trajectory of the light ray, that is in t* and 
k. Therefore, the function A(r, tq) must be viewed as depending on four variables A = A(s, sq, t*, fe). Accounting for 
these remarks the derivative along the light ray reads as follows 
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where the unit vector K is defined in (6.30) and where we explicitly show the dependence of function A(t, tq) on all 
parameters which implicitly depend on time. 

The time derivative of vector k is calculated using the approximation k = K and formula (6.30) where the 
coordinates of the source of light, a;o(io)i and of the observer, x{t), are functions of time. These derivatives are 

dk' _ (fc X (y X k))' 9fc' _ (fc X (Vo X k))' 

dt R ' dto R ' ^ ^ 

where R = |a; — a;o| is the distance between the observer and the source of light. 

Time derivatives of retarded times s and sq with respect to t and Iq are calculated from their definitions, s = t — r 
and So = ~ ^Oj where we have to take into account that the spatial position of the point of observation is connected 
to the point of emission of light by the unperturbed trajectory of light, x ^ Xq + k {t — t^) . More explicitly, we use 
for the calculations the following relationships 

s = t - \xo{ta) + k {t - to)\ , (7.38) 
so = h-\xo{tQ)\ , (7.39) 

where the unit vector k = fe(t, to) must be also considered as a function of two arguments to with its derivatives given 
by equation (7.37). We notice that relationship (7.38) combines two null cones related to propagation of gravitational 
field from the isolated system and to that of an electromagnetic wave from the source of light to observer. Equation 
(7.39) is just a null cone corresponding to propagation of the gravitational field from the isolated system to the point 
of emission of light. Calculation of infinitesimal variations of equations (7.38), (7.39) immediately gives a set of partial 
derivatives 

= 1 - fc • AT - (fc X y) • (fc X AT) , (7.40) 



dt 
ds 

dsp 
dto 

dt 



il-k-Vo)ik-N) , (7.41) 
1 - Fo • ATo , (7.42) 
, (7.43) 



where — x^/r and Nq — Xq/tq are the unit vectors directed from the isolated system to observer and to the source 
of light respectively. . 

Time derivatives of the parameter t* = t^ — k ■ a;o(io)j where again fc = fc(t, to), read 

The last terms of order \^\/R in both formulas are generated by the time derivatives of the vector fc. 

In what follows we shall restrict ourselves with the case of static observer and the source of light, that is velocities 
V = Vq = 0. Taking into account this restriction in equations (7.37)-(7.44) and leaving only linear with respect to 
the universal gravitational constant G terms, allow us to simplify equation (7.36) 

Function A = A(s, sq, t*, fc) is defined by the integral (7.2) which (in the case of static observer and the source of light) 
depends on the retarded times s, so and the time of the closest approach t* through the only arguments y = s — t*, 
2/0 = So — t* and t* , that is 

A(s, so, t*, fc) = A(y, yo, t*) = ^{y, t*) ~ ^(yo, t*) , (7.46) 

where the relativistic perturbation of the eikonal, tp = — fc'S, is defined in accordance with equation (7.6). Therefore, 
the partial derivatives with respect to s and t* from the argument y, and those with respect to sq and t* from yo in 
equation (7.45) are nullified and it is reduced to simpler form 

dto ^ d^P{y,t*) d4,{yo,t*) , . 

dt dt* dt* ' ^ ' ^ 
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The partial time derivative from ip(y,t*) is found by differentiation of relationships (6.2)-(6.4) 
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When deriving these equations we assumed that the time evolution of both the mass M and the angular momentum 
5* of the isolated system can be neglected. In opposite case, the right side of equations (7.50) and (7.51) would 
contain also time derivatives from the mass and the angular momentum. Partial derivative from function ^j{yo, t*) is 
obtained from equations (7.48)-(7.51) after replacements y — > yo, t tg, and r tq. 

Provided that the frequency shift z is given by equation (7.31) we find that in the case when both observer and 
the source of light are static with respect to the ADM-harmonic reference frame, the frequency (energy) of a photon 
propagating through gravitational potential changes in accordance with (exact) relationship 
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In the linear with respect to G approximation equation (7.52) yields 
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(7.53) 



where Siy = f — Vq. One notices that the time of the closest approach t* enters equation (7.53) explicitly. At the first 
glance the reader may think that it plays a significant role in calculations and must be known in order to calculate the 
gravitational shift of frequency. However, the partial derivative with respect to t* must be understood in the sense of 
equation (4.27) which makes it evident that the time derivative with respect to t* taken on the light-ray path is, in 
fact, the time derivative with respect to time t before using the light-ray substitution. Thus, the gravitational shift 
of frequency can be recast to the next practical form 



5v_ 



hQo{t,x) - hoo{to,xo) 
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dt 



dto 



(7.54) 



where ^/'(t, a?) and ip{tQ,XQ) are relativistic perturbations of the electromagnetic eikonal taken at the points of ob- 
servation and emission of light respectively. These time derivatives are, of course, the same as shown in equations 
(7.48)-(7.51). Equation (7.54) has been derived by making use of the definition (7.31). It is straightforward to prove 
that the definition (7.32) brings about the same result. Indeed, in the case of the static observer and the source of light 
their 4- velocities are simply u" = {dt/dT, 0, 0, 0) and Uq = {dto/dTo, 0, 0, 0), and the wave vector = ^{ka + daip) (see 
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equation (4.33)) where lo = 27ri/oo is a constant frequency which is conserved along hght ray because of the equation 
of the paraUel transport (4.13). Substituting these definitions to equation (7.32) leads immediately to equation (7.54) 
as expected 

Physical interpretation of the relativistic frequency shift given by equation (7.54) is straightforward although the 
calculation of different components of the time derivatives from the eikonal are tedious. The term in the square 
brackets in the left side of equation (7.54) reads 



hoo{t,x) - hao{tQ,XQ] 



GM _ GM 
r tq 



(7.55) 



and represents the difference between the values of the spherically-symmetric part of the Newtonian gravitational 
potential of the isolated system taken at the point of observation and emission of light. The time derivatives of the 
eikonal depending on the mass and spin multipole moments of the isolated system are given in equations (7.50), (7.51). 
Scrutiny examination of these equations reveal that these components of the frequency shift depend on the first and 
higher order time derivatives of the multipole moments and would vanish in the stationary case. The gravitational 
frequency shift contains the gauge-dependent contribution as well. This contribution is given by equation (7.49) and 
can be calculated by making use of equations from appendix C. The result is as follows 
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Time derivatives dt^{k^w^ — w°) and dt'{k'^(p^ — (f'^) can be obtained from equations (7.56), (7.57) after making 
replacements t to, r ^ vq, and r ^ tq. 



7.4. Gravity-Induced Rotation of the Plane of Polarization of Light 

Any stationary or time-dependent axisymmetric gravitational field in general induces a relativistic effect of the 
rotation of the polarization plane of an electromagnetic wave propagating in this field. This effect is similar to 
Faraday's effect in electrodynamics [126]. The Faraday effect is caused by the presence of magnetic field while the 
gravity-induced rotation of the plane of polarization of light is caused by the presence of the, so-called, gravito- 
magnetic field [127]. This gravito-magnetic effect was first discussed by Skrotskii ([108, 109] and shortly afterwards 
by many other researches (see, for example, [110, 111, 128-132] and references therein). Recently, we have studied 
the Skrotskii effect caused by a spinning body moving arbitrarily fast and derived Lorentz-invariant expression for 
this outstanding effect [133] in the linearized approximation of general relativity containing all terms of order v/c 
coupled with the universal gravitational constant G. In the present paper we further generalize our theoretical results 
along with those of previous authors regarding the case of an isolated system emitting gravitational waves of arbitrary 
multipolarity. 
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In order to simplify description of the Skrotskii effect we shall consider the parallel transport of the reference tetrad 



defined by equations (4.55) at each point along the light ray. We shall assume for simplicity that at the past null 
infinity the tetrad has the property that e'jjj(— (X))9°j-) = in accordance with equation (4.50). The parallel transport 
of this tetrad is defined by equation (4.66) which makes clear that in an arbitrary point on the light-ray trajectory 
e'^^^-^{t,x) = 0{hap). We are interested only in solving equation (4.66) for the vectors e"^^ {n = 1,2) that are used in 

the description of the polarization of light. The propagation equation for the spatial components e*^^-) of these vectors 
are directly obtained from equation (4.66) which yields 
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Here we have defined the angular velocity of gravito-magnetic drag by fi* defined as 



(7.59) 



This quantity describes the rate of the rotation of the plane of polarization of electromagnetic wave caused by the 
presence of the gravito-magnetic field. As soon as equation (7.58) is solved the time component e^^^^ is obtained from 
the orthogonality condition, lae"(n) — 0; which implies that 
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where the relativistic perturbation of the light-ray trajectory is given in equation (6.15). 

Let us decompose vector into three components that are parallel and perpendicular to the unit vector fc*. We 
can use in the first approximation the well-known decomposition of the Kroneker symbol in the orthonormal basis 



(7.61) 



where (a, b, k) are three orthogonal unit vectors of the reference tetrad at infinity. Decomposition of fi* is, then, given 

by 



= {a- n)a^ + (b • 0)6' + (fe • n)k^ 



(7.62) 



Taking into account that vectors e^j-j^^ = + 0{hap), e*^2) = + 0{hap) one obtains the equations of the parallel 
transport of these vectors in the next, more convenient for integration, form 
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where equalities Sijie'^i^-^^k^ = —&^(2) + 0{haf}) and Sijie^f^^)^^ — + Oihaji) have been used. 

Solutions of equations (7.63), (7.64) in the linear approximation with respect to the angular velocity fJ* reads 
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where the second terms in the right sides of equations (7.65), (7.66) preserve orthogonality of vectors e'^j^j and 

in the presence of gravitational field while the last terms describe a small rotation (the Skrotskii effect [108, 109]) of 

each of the vectors at the angle 
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about the direction of light propagation, fc, in the local plane of vectors e(i) and 6(2). It is worth noting that the 
Euclidean dot product fc • fl can be expressed in terms of partial differentiation with respect to the vector of the 
impact parameter. This can be done by making use of equation (4.28) and noting that Eijpk^k^ = 0, so that 

k Vt^ \k"yejf,qdqKp , (7.68) 

where the hat over spatial indices denotes the projection onto the plane orthogonal to the propagation of light ray, 
for instance, A* = P^A^ ■ Hence, the transport equation for the angle cfi assumes the following form 

g = h^k^e^pgdgh^p . (7.69) 

This equation can be split in three linearly-independent parts corresponding to contributions of the gauge, (/'(o, the 
mass, (^(Af), and the spin, (f>(s), modes of the gravitational multipolar field to the overall value of the Skrotskii effect. 
Specifically, we have 

= 0(G) + 0(M) + 0(S) + 00 , (7.70) 

where (f>o is constant angle characterizing initial orientation of the polarization ellipse in the plane formed by the e^^-) 
and 6(2) vectors. 

The gauge-dependent part of the Skrotskii effect is easily integrated so that wc obtain 

Aa, = lk^e„,d,{wP+xn , (7.71) 

where the gauge functions = wf^j +u'fs) and x* = xf*/) +X(s) are given in equations (3.20) and (C.IO) respectively. 
More specific calculations reveal that the gauge-dependent terms enter the right side of equation (7.71) in the following 
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form 
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The reader can notice the integrals in the right side of equations (7.72) and (7.73). The integrals are actually should 
not be calculated explicitly since there is sufficient number of partial derivatives in front of them which makes the 
integration unnecessary in correspondence with the rules of differentiation of such integrals worked out in section 5. 
The Skrotskii effect due to the mass multipoles of the isolated system is given by 
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The gravitational field of the spin multipoles rotates the polarization plane of the electromagnetic wave at the following 
angle 
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Integrals in equations (7.76) and (7.77) are eliminated after taking one partial derivative so that we do not need to 
integrate. 
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8. OBSERVABLE RELATIVISTIC EFFECTS IN THE GRAVITATIONAL LENS APPROXIMATION 

Here we will consider a limiting case of gravitational lens approximation when the impact parameter d of the light 
ray with respect to an isolated system is much smaller than distance rg from the isolated system to the source of light 
and distance r from the isolated system to observer (sec Figs. 3). 



In the case of a small impact parameter of the light ray with respect to the isolated system the near-zone gravitational 
field of the system strongly affects propagation of the light ray when the light particle (photon) moves in the vicinity of 
the system so that the effects due to the gravitational waves are suppressed [82, 83]. In what follows, we assume that 
the impact parameter d of the light ray is small comparatively with both distances r and rp, that is, d <^ min[r, rg] 
(sec Fig. 3). This assumption allows us to introduce two small parameters: e = d/r and = d/vQ. If the light ray 
is propagated through the near-zone of the isolated system one more small parameter can be introduced, specifically, 
the parameter ex = d/X, where A is the characteristic wavelength of the gravitational radiation emitted by the system 
proportional to the product of the speed of gravity and the characteristic period of oscillations of the isolated system. 
One emphasize that the parameters e, Eq and e\ are not physically correlated. Parameter e\ is used for the post- 
Newtonian expansion of the observed effects - deflection of light, time delay, etc. This expansion can be also viewed 
as a Taylor expansion with respect to the parameter u/c, where v is the characteristic speed of motion of matter 
composing the isolated system and c is the speed of propagation of gravity (equal to the speed of light . If the light 
ray does not enter the near zone of the system the parameter Ex is not small and the post-Newtonian expansion can 
not be performed. 

Small-impact-parameter expansions for variables y and ya yield 



8.1. Useful Asymptotic Expansions 




(8.1) 



(8.2) 



and 




(8.3) 




(8.4) 



where the numerical coefficient of the expansions is 




(8.5) 



Time variables t and io are expanded as follows 




(8.6) 



(8.7) 



where Ck is given by equation (8.5) and t* is the time of the closest approach of light ray to the barycenter of the 
isolated system. 
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Using (8.6) and (8.7) we can write down the post-Newtonian expansions for functions of the retarded time t — r as 
follows 



Fit-r) = J2 



k=0 



(-l)'^d'^ ( e 



kl 



.2k-l 



k=2 



-2fc-l 



fc! \ 2 

fe=0 \ fc=2 



F^''\t* - 2ro) = F{t* - 2ro) + eo-F{t* - 2ro) + O {elel) . (8.9) 



Wc notice that convergence of the post-Newtonian time series given above depends, in fact, not just on one parameter 
ex which is typical in the post-Newtonian celestial mechanics of extended bodies, but on the product of two parameters. 
Therefore, the convergence requires satisfaction of two conditions 



EoEa < 1 , 



(8.10) 
(8.11) 



where the numerical value of the parameter must be taken for the smallest wavelength, Amin, in the spectrum 
of the gravitational radiation emitted by the isolated system. Conditions (8.10), (8.11) ensure convergence of the 
post-Newtonian series in (8.8) and (8.9) respectively. If the source of light rays and observer are at infinite distances 
from the isolated system then e ~ and eo — 0, and the requirements (8.10), (8.11) are satisfied automatically, 
irrespective of the structure of the Fourier spectrum (5.5) of the multipole moments of the isolated system. In real 
astronomical practice such an assumption may not be always satisfied. In such cases it is more natural to avoid the 
post-Newtonian expansions of the metric tensor and/or observable effects. 

If we assume that the mass and angular momentum of the isolated system are conserved the asymptotic expansions 
of integrals (5.1), (5.2) of the stationary part of the metric tensor have the following form 
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(8.12) 
(8.13) 
(8.14) 
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Some estimates for the derivatives with respect to ^' and r of some functions: 
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where to > 2, n > 0. Here and in the subsequent expressions (8.17) - (8.23) m and n depend on I] we give the estimates 
from below for m and n for ^ > 1. 
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where to > 2, n > 0; 
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(8.18) 
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where to > 2, n > 0: 
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where to > 1, n > 0; 
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where to > 2, n > 0: 
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where to > 3 and ri > 0. 

Having written down the auxihary asymptotic expansions we can consider the observable relativistic effects in the 
gravitational lens approximation. 

8.2. Asymptotic expressions for observable effects 

Here we write out the expressions for relativistic effects neglecting all the first order terms with respect to e. 
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The observable unit vector in the direction " observer - source of light" is given by the expression 



(8.26) 
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where we have neglected the quantities (tq,^) and 7* (t, ^). The vector {t,^), characterizing the deflection of 
hght, and the relativistic correction /?' (t, ^) are given by the expressions 

a'(T,0-«j„,(r,4) + a;„(r,0, (8.28) 

=4X4lnrf+ (8.29) 



where 
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where 
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/3^(T,0-/3(M,(r,0+/3:.,(r,0, (8.31) 

/3:m,(^,«)=-^«U +«(£«), (8.32) 

/3^,)(T,|)=-^<„+0(e«). (8.33) 

We note that using the expressions (8.8) and (8.9), we can rewrite (8.25), (8.26), (8.29) and (8.30) expressing 
functions at the moment of time (t — r) in terms of their values at the moment of the closest approach of the photon 

and the gravitating system t* . This is accomplished simply by substituting the value t* instead of {t — r) in the 
expressions (8.25), (8.26), (8.29) and (8.30) because the corrections will be of the first order with respect to the 
parameter e (see Eq. (8.8)). 

A (r. To) ^2M ln(4rro) - 4^/^, (8.34) 

a*(T,0=4a,VA (8.35) 
where '0 - is the gravitational lens potential, which has the form 

00 1-2 ^_^y^p 
1=2 p=0 

df.lA, it*)k^ai...apdaj,+ i...a,> + 

^jbakjShda — 

00 ;-i 



1=2 p=0 

^?'SbA^_^it*)k<a^...ap^ap+^...a,_^> ^lnd + 0{eex). 

This expression takes into account multipole moments of all orders in the geostationary case and presents a general- 
ization of the corresponding results obtained in [7] , [84] and [83] . 

The expressions for the angle of rotation of the polarization plane of light in gravitational lens approximation have 
the form 

(r ) - 0<.„ (to) = - 4 5] ^ Ci il-p,p)j^x 

i=2 p=0 

X 4^*ej6a,2^byii_i(^*)^«ii-ap^ap+i-a!>j lnrf + 0(eeA), (8.37) 
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9. OBSERVABLE RELATIVISTIC EFFECTS IN THE PLANE GRAVITATIONAL-WAVE 

APPROXIMATION 



9.1. Plane- wave Asymptotic Expansions 

In the plane gravitational wave approximation we assume that A ^ min[r,ro] (wave zone condition) and 
R <^ min[r, ro] (plane wave condition) (Fig. 2). We introduce small parameters S\~X/r, S — R/r and Sq = R/r^ 
, where A is the characteristic wavelength of the gravitational radiation emitted by the isolated astronomical system. 
The relationship between r, rg and R can be written as follows 



rl -2rRcose + R^ = r^{l - 25 cosO + 5^) , 



(9.1) 



where 9 - is the angle between the directions "observer - the source of light" and "observer - the gravitating sys- 
tem" (Fig .2). From the previous expression it follows that 



ro = r(l - 6 cose + 0{5'^)) , 

2\ 



— = -(1 + 6 cose + 0(5')) 
ro r 

For the variables t and tq we have 

T = r cos e, To = T — R ^ r cos e — R. 
Quantities d and y can be expressed as follows 

d ~ r sin 9 , 

y = T — r = — r(l — cos9) . 
The instants of time {t — r) and {to — ro) are related to each other as follows 

to — ro = t — r — R{1 — cos 9) . 
Using this expression we can write the Taylor expansion for the functions of retarded time s — t — r 

F{to - ro) = F{t - r) - R{1 - cos9)F{t - r) + O {(R/X)^) • 



(9.2) 
(9.3) 

(9.4) 
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Let us now write out the asymptotic expressions for the derivatives with respect to and t of some functions 
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where = x'/r. Expressions (9.13) can be proved by induction. Expressions for functions taken at the retarded 
instant of time sq = to — tq can be obtained from equations (9.13) by substitution to, tq and 9q instead of t, r and 
respectively. 

Finally, we notice that the impact vector can be decomposed as follows 



f = r (TV* - fc'cose) . 
This form of the impact vector ^' is useful in subsequent approximations. 



(9.14) 



9.2. Asymptotic Expressions for Observable Effects 



In this section we give expressions for the relativistic effects of the time delay and bending of light in the gravitational 
plane- wave approximation. In this approximation we neglect all terms of order 5^, Sq and (5^. For the time delay one 
has 



A - A (t,to) + A (t,to), 

(M) (S) 



where 



■E 



(jvf) r 1 — cos {) ^-^ i! 

1=2 



-, TT 



iijAi^^ih - ro) 



TT 



A r, To = - 2- 

(s) 1 — cosw \r ro 



(9.15) 

(9.16) 
(9.17) 



4fc,: 



i-iyi 



. (l + l 



1- cos {i + iy. 



TT 



H TT 



^0 



- ,aAi 



Here the transverse - traceless (TT) part of the tensors depending on the multipole moments is taken with respect to 
the direction N"^ . Taking into account expression (3.8) for the components of the metric tensor hij we can re- write 
expressions (9.15) - (9.17) for the time delay as follows 



A = A (r, To) + A (r, tq) 

[M) (S) 



(9.18) 



r 1 — cos 



In a similar way we obtain expressions for the observable direction from observer to the source of light perturbed by 
the gravitational wave 



where 



2 1 — cos 



[{cose- 2)F + N'] hl^{t - r) + k''hjy^{t - r), 



(9.19) 

(9.20) 
(9.21) 
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In expression (9.19) we have neglected the quantities /3* (t, ^). Expressions (9.18) - (9.20) were obtained in [83] in 
the spin-dipole mass-quadrupole approximation only. This paper has generalized this result for the case when hjj^ 
contains all mass and spin multipole moments of arbitrary orders. 

In the case when the distance between observer and the source of light is much smaller than the wavelength of the 
gravitational waves i? <C A, expression (9.18) is reduced to a well known result [1, 2] for gravitational wave detectors 
located in a wave-zone of an isolated system 

\h,hjnt - r) + 0{5^) + o(^\ , (9.22) 



where AR = c A (r, tq). 

Relativistic rotation of the polarization plane of light (the Skrotskii effect) in the gravitational plane-wave approx- 
imation assumes the next form 



ir) - Hro) ^ -,'-Y^hl-it - .) - -, \_^\ hj-it, - ro) . (9.23) 



APPENDIX A: APPROXIMATE EXPRESSIONS FOR THE CHRISTOFFEL SYMBOLS 

Making use of the general definition of the Christoffel symbols given in equation (4.14) and applying the expansion of 
the metric tensor (3.1) result in the approximate linearized post-Minkowskian expressions for the Christoffel symbols 

= -ldthoo{t,x) , 

r°- = -^[d^hoj{t,x) + ^Jho^{t,x) ^ dth,j{t,x)] , (A.3) 

Too - dthQ,{t,x)~^d,hoo{t,x) , (A.4) 

Toj = ^[djhoi{t,x) - dihoj{t,x) + dthij{t,x)] , (A. 5) 

Tjp = ^[djhip{t,x) + dphij{t,x) - dihjp{t,x)] . (A. 6) 
These expressions are used for the calculation of the right-hand side of the null geodesic equation (4.13). 

APPENDIX B: THE METRIC TENSOR IN TERMS OF THE LIGHT-RAY PARAMETERS | AND r 

Expressions for perturbations of the metric tensor (3.5) - (3.8) in terms of the parameters and r have the 
following form follows 
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where 
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All quantities in the right side of expressions (B.5)-(B.9), which are explicitly shown as functions of x\ r = and t, 
must be understood as taken on the unperturbed hght-ray trajectory and expressed in terms of = |^|, t and t* in 
accordance with equations (4.20), (4.22). For example, the ratio X^~'-'(f — r)/r in equation (B.5) must be understood 



as 



47"^ - _ I a'" + v^^Trf^) 

and the same replacement rule is applied to the other equations. 
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APPENDIX C: THE GAUGE FUNCTIONS 

Gauge functions w", generating the coordinate transformation from the canonical harmonic coordinate system to 
the ADM-harmonic one arc given by equations (3.19), (3.20). They transform the metric tensor as follows 



(C.l) 



where /i^™' is the canonical form of the metric tensor in harmonic coordinates given by equations (3.5)-(3.8) and hap 
is the metric tensor given in the ADM-harmonic coordinates by equations (3.21)-(3.24). 

The gauge functions taken on the light-ray trajectory and expressed in terms of the variables ^ and r can be written 
down in the next form 
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where h qq''(u, t, ^), hX'^ {u,t, ^) and hX'^{u,T,^) are defined by the Eqs. (B.5), (B.6) and (B.7) after making use of 

(J\/) (M) (S) 
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the substitution t* — > u. It is worth noting the following relationships 
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which are helpful in calculation of the gravitational frequency shift. 

Linear combination k'^ip'^ — ip^ of the gauge-dependent functions introduced in equation (4.41) is given by the 
expressions 
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Gauge-dependent term generated by equation (7.69) for the rotation of the plane of polarization of electromagnetic 
wave is a pure spatial vector that can be decomposed in two linear parts corresponding to the mass and spin 
multipoles: 
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FIG. 1: Penrose's diagram shows a worldline of the isolated system (binary star) originating at past timelike infinity and 
ending at the future timelike infinity 7^. Both light and gravity propagate along null geodesies going from the past null infinity 
J~ to the future null infinity . A particular light geodesic is emanating from the event (to, xq) and ending at the event (t, x). 
Extrapolation of this light geodesic to the past null infinity defines the null wave vector fc" = (fe°,fe) of the electromagnetic 
wave under consideration. 
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To observer 



Vector field of light rays 



To observer 



Plane of the sky 



FIG. 2: Astronomical coordinate system used for calculations. The origin of the coordinate system is at the center-of-mass of 
the source of gravitational waves. The bundle of light rays is defined by the vector field fc*. The vector = — fc* + 0{c~^) 
is directed from observer towards the source of light. The vector Kq is directed from the observer towards the source of 
gravitational waves. We use in the paper the equalities Kq = —A*'' = —x''/r, where are the coordinates of the observer with 
respect to the source of gravitational waves, and r = \x\. The plane of the sky to the vector is not shown. 



60 




Source of gravitational waves 



FIG. 3: Relative configuration of observer (O), source of light (S), and a localized source of gravitational waves (D). The source 
of gravitational waves deflects light rays which are emitted at the moment to at the point S and received at the moment t at 
the point O. The point E on the line OS corresponds to the moment of the closest approach of light ray to the deflector D. 
Distances are OS = R, DO = r, DS — ro, the impact parameter DE = d, OE = r > 0, ES — tq = t — R < 0. The impact 
parameter d is small in comparison to all other distances. 
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FIG. 4: Relative configuration of observer (O), source of ligfit (S), and a localized source of gravitational waves (D). The source 
of gravitational waves deflects light rays which are emitted at the moment to at the point S and received at the moment t at 
the point O. The point E on the line OS corresponds to the moment of the closest approach of light ray to the deflector D. 
Distances are 05* = R, DO — r, DS = ro, the impact parameter DE — d, OE = r — rcosO, ES = to = t — R. The distance 
R is much smaller than both r and ro. The impact parameter d is, in general, not small in comparison to all other distances. 



